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Abstract 

Sound  propagation  and  the  Initial  value  problem  in 
kinetic  theory  are  Investigated.   Linearized  forms  of  certain 
kinetic  models  suggested  by  Gross  and  Jackson  are  employed. 
The  most  general  model  contains  three  relaxation  times  and 
Is  capable  of  producing  Euler,  Navler-Stokes ,  Burnett  and 
thirteen  moments  equations  for  smooth  phenomena.   A  study 
of  dispersion  relations  is  made  and  some  novel  features  are 
uncovered.   One  finds  that  kinetic  models  are  unable  to 
describe  phenomena  of  higher  than  a  certain  wave  number,  the 
latter  depending  on  the  model  chosen.   It  compensates  for  this 
by  introducing  an  interesting  but  unphysical  dispersion 
picture  for  higher  wave  numbers.   It  is  further  suggested  by 
one  of  the  models  that  the  phase  speed  of  sound  waves  achieves 
a  maximum  value. 

Existence,  uniqueness,  and  boundedness  of  solutions  of 
the  initial  value  problem  are  shown  for  any  model  equations. 
It  is  further  shown  that  asymptotically  the  solutions  become 
hydrodynamical . 
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Sound  Propagation  According  to  Kinetic  Models 


Introduction. 

Sound  propagation  by  its  very  nature  is  a  molecular 
effect.   Yet  it  was  not  until  the  work  of  Wang  Chang  and 
Uhlenbeck   that  sound  propagation  was  studied  from  a  molecu- 
lar viewpoint.   Prior  analytical  attempts  were  essentially 
hydrodynamical  and  its  extensions.   By  the  latter  we  shall 
mean  the  equations  obtained  by  means  of  the  Chapman-Enskog 
procedure.   In  a  sense,  the  work  of  Wang  Chang  and  Uhlenbeck 
is  also  macroscopic.   They,  in  fact,  considered  their  method 
a  Chapman-Enskog  process  and  labelled  their  calculations  as 
Euler,  Navier-Stokes ,  Burnett,  etc.   These  appellations  can 

be  more  appropriately  given  in  terms  of  the  moments  equa- 

2  "5 
tions  of  H.  Grad.  '    This  was  pointed  out  in  a  recent 

study  by  one  of  the  authors.    From  the  results  of  the  latter 

report  it  is  clear  that  a  complete  dispersion  theory  of  sound 

propagation  using  the  Boltzmann  equation  is  made  unfeasible 

by  the  difficulty  in  obtaining  results  for  large  wave  number 

phenomena.   It  is  our  goal  to  obtain  a  description  valid 

throughout  the  wave  number  range. 

To  accomplish  this  the  Boltzmann  equation  is  abandoned 

in  favor  of  certain  kinetic  models  suggested  by  Gross  and 

Jackson.-^  We  will  neither  obtain  these  equations  from  the 

Boltzmann  equation  nor  relate  it  to  that  equation.   The  former 
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discussion  may  be  found  in  reference  5  and  the  latter  in 
reference  4.   We  will  be  content  with  knowing  that  our  most 
general  equation,  given  in  Section  2,    shares  many  common 
features  with  the  Boltzmann  equation.    Although  we  do  not 
relate  the  model  equations  to  the  Boltzmann  equation,  we  can 
look  in  the  other  direction  and  relate  them  to  macroscopic 
equations.   We  will  find  that  our  general  equation,  which 
contains  three  constants  (and  is  referred  to  as  the  three 
relaxation  time  model)  is  capable  of  producing  the  Euler, 
Navier-Stokes,  Burnett,  and  thirteen  moments  theories.   In 
the  course  of  our  investigations  we  will  introduce  several 
approximate  equations  which  will  be  referred  to  as  isosteric 
(mass  preserving),  isothermal  (constant  temperature),  and  the 
single  relaxation-^  cases.   In  no  case  will  we  have  the  latitude 
of  our  general  equation;  in  fact  the  latter  approximation, 
which  is  the  most  general,  is  only  capable  of  furnishing  the 
Euler  theory.   It,  therefore,  seems  necessary  to  choose  a  case 
as  general  as  the  three  relaxation  model  if  more  than  qualita- 
tive agreement  with  physics  is  sought. 

In  their  well-known  paper,  Bhatnager,  Gross  and  Krook  , 
among  other  things,  consider  sound  propagation.   Their  equa- 
tion, which  is  now  referred  to  as  the  single  relaxation  model 
mentioned  above,  is  considered.   Actually,  they  only  carry 
out  the  analysis  for  what  we  have  called  the  isothermal  case. 
Our  methods  and  results  are  sufficiently  different  to  warrant 
a  reconsideration  of  this  case.   The  most  notable  difference 
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is  the  exhibition  of  a  part  of  the  spectrum  which  was  not 
previously  discussed.   It  can  be  legitimately  referred  to  as 
the  analytic  continuation  of  the  dispersion  relation  but, 
nevertheless,  is  still  part  of  the  spectrum.   During  the  •. 
course  of  the  discussion  an  argument  is  given  which  shows  -    ■■■ 
that  this  analytical  continuation  of  the  spectrum  has  no 
physical  basis.   This  spurious  part  of  the  spectrum  arises 
because  a  model  equation  regards,  in  a  sense,  phenomena 
past  a  certain  wave  number  as  free-flow.   The  dispersion 
relation  of  a  model  equation,  therefore,  may  only  be  used 
with  confidence  for  phenomena  of  low  enough  wave  number.   The 
threshold  value  of  wave  number  depends  on  the  model  chosen. 

In  spite  of  the  approximate  nature  of  the  model  equations, 
actual  results  from  the  dispersion  relation  are  difficult  to 
obtain  except  for  small  wave  numbers.   However,  the  disper- 
sion relation  does  appear  in  closed  form.   Although  the 
dispersion  relations  are  amenable  to  numerical  calculation, 
only  a  few  calculations  were  made.   A  report  now  in  progress 
will  contain  an  extensive  numerical  study  of  the  dispersion 
relations.   The  few  calculations  do  show  the  possibility  of  a 
maximum  phase  speed  --  it  is  implied  by  the  single  relaxation 
model.   This  will  not  be  explored  but  must  await  the  future 
numerical  work. 

The  entire  report  is  attuned  to  solving  an  initial  value 
problem  in  the  absence  of  boundaries.   The  boundary  value 
problem  is  physically  more  desirable  but  unfortunately  more 
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difficult  to  deal  with.   Inasmuch  as  the  dispersion  relation 
Is  concerned,  the  boundary  value  problem  requires  solving 
for  wave  number  In  terms  of  the  frequency,  whereas  the  pure 
Initial  value  problem  requires  the  Inverse  relation.   In  the 
forthcoming  numerical  work  the  former  case  will  also  be  con- 
sidered. 

Since  the  physical  basis  of  kinetic  model  equations  is 
somewhat  obscure,  a  mathematical  exploration  of  linearized 
model  equations  is  made  in  the  last  section.   The  model  equa- 
tions like  the  Boltzmann  equation,  are  integro-differential 
equations  and  lend  themselves  very  nicely  to  the  standard 
procedures  used  in  showing  the  existence  and  uniqueness  of 
solutions.   The  latter  properties  are  shown  for  all  time 
under  rather  Innocuous  restrictions  on  the  Initial  data.   It 
is  in  fact  shown  that  the  solutions  are  exponentially  bounded 
in  time,  which  validates  the  use  of  Laplace  transforms. 
Further,  using  the  results  of  the  dispersion  theory,  one  can 
show  that  the  solutions  are  uniformly  bounded  in  time  and 
behave  asymptotically  as  hydrodynamics. 

Acknowledgement:   During  the  course  of  working  on  this  material, 
both  authors  profited  from  discussions  with  many  colleagues. 
Our  interpretation  of  results  and  our  point  of  view  has  been 
notably  enriched  by  discussions  with  J.  Berkowitz,  C.  Gardner, 
H.  Grad,  and  H.  Weitzner. 
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1 .   Hydrodynamic  Preliminaries.  \ 

Before  developing  the  kinetic  theory  description.  It  Is 
advantageous  to  first  consider  the  hydrodynamic  description 
of  sound  propagation.   This  will  be  helpful  for  later  com- 
parison.  A  convenient  starting  point  Is  the  linearized 
one -dimensional  Navler-Stokes  equations, 


(^•1)       s^ 


+  |H=  0 


(-]■->)         aT     2  su  _  .  a^    ^ 


All  quantities  (which  are  perturbations  from  equilibrium) 
occurlng  In  (1.1-3)  and  the  temporal  coordinate   t,   and  the 
spatial  coordinate  x  have  been  made  dlmenslonless  with 
respect  to  the  mean-free-path  L  and  a  collision  time 
defined  by 


(lA) 


m^ 


R  Is  the  gas  constant  and  T   the  unperturbed  temperature 
Denoting  all  the  unperturbed  quantities  by  zero  subscripts, 
and  their  dimensional  perturbations  by  circumflexes,  the 
normalization  Is  given  by 
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(1.5)     P  -  p/pQ  ,   density 


(1.6)     p  =  P/Po  '    pressure 


(1.7)     T  =  T/T   ,    temperature 


(1.8)      u  =  u//RT^  ,  velocity, 


Space  has  been  normalized  by  L  and  time  by   t.   The  two 
dimensionless  quantities   e   and   6   are  given  by 


(1.9)      e  = 


4v 


3l/rt^ 


(1.10)     5  =  — ^- 

l/rt 


where   v   is  the  kinematic  viscosity  and   k  the  thermo- 
metric  conductivity.   Further,  the  ratio  of  specific  heats 
(y  =  c  /c  )   has  been  taken  as  5/3-   Under  the  linearization 
and  normalization  the  "gas  law"  is  given  by, 

(1.11)    p  =  p  +  T  . 

Equation  (l.ll)  will  be  used  to  eliminate  the  pressure 
in  the  sequel. 

It  will  also  be  of  interest  for  us  to  now  consider  the 
physically  unreasonable  but  simpler  case  of  isothermal 
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propagation  gotten  by  taking  T  =  0, 


(1.12) 


(1.15) 


dt 

0 

dx 

e 

This  case  may  be  thought  of  as  coming  from  equations  (1.11-5) 
by  making  the  heat  conductivity  unbounded. 

We  now  consider  the  Initial  value  problem  for  (1.12-13) 
In  the  absense  of  boundaries.   This  will  be  the  only  problem 
considered  j.n  this  report  arid  we  take  this  Initial  value  problem 
and  sound  propagation  as  being  synonymous.   On  taking  a 
Laplace  transform  In  time  and  a  Fourier  transform  In  space 
and  using  matrix  notation  the  Isothermal  system  Is 


(1.14) 


-Ik 


■Ik 


a  +  ek 


Initial 
Conditions 


/ 


a   and   k  denote  the  Laplace  and  Fourier  variables  respec- 
tively.  Since  we  will  not  be  Interested  In  solving  any 
particular  Initial  value  problem  the  right  hand  side  of 
equation  (l.l4)  Is  left  unspecified.   We  are  Interested  In 
the  roots  of  the  determinant  of  the  matrix  of  (l.l4).  I.e., 
In  the  dispersion  relation.   Strictly  speaking,  therefore,  we 
shall  be  Investigating  the  class  of  possible  plane  wave  solu- 
tions . 
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The  roots  of  the  isothermal  dispersion  relation  are  easily 
computed  and  given  by 

(1.15)    a^  =  2 

2 
If  this  Is  plotted  in  the  complex  a-plane  with  k   as  a 

parameter,  we  get  the  sketch  of  Figure  1.   There  are  two 

o propagating  branches  starting  at  the  origin,  one  to  the  right 

and  the  other  to  the  left.   For  a  sufficiently  large  value 

of   k  (k  =  4/e  )   the  two  branches  return  to  the  negative 

real  axis,  as  is  indicated  in  Figure  1,  and  then  continue 

as  two  non-propagating  diffusion  modes.   One  stops  at  a  finite 

point  while  the  other  continues  to  negative  Infinity.   The 

details  of  the  latter  effect  do  not  really  concern  us  since 

from  the  point  of  view  of  kinetic  theory,  hydrodynamical 

equations  can  only  be  considered  for  small  wave  number  k. 

In  this  limit 


ek^    ie^k^  ,  ^,,4 


(1.16)  a^  =  Ik  -  ^  -  ^■V^+  O(k-) 

2     2  3 

(1.17)  a_  =  -ik —  +  — g —  +  0(k  ) 


As  will  always  be  the  case,  the  roots  occur  in  conjugate  pairs 
Splitting  a  ,  for  instance,  into  real  and  Imaginary  parts. 


(1.18)    o^  =  al  +   lal    , 
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we  e:et 


(1.19)     a^=  -S|!+  oCk") 


(1.20)     !±  =  Vp  =  1  -  4^  +  0(k'*) 


The  quantity  (I.19)  gives  the  attenuation  rate  of  a  plane 
wave  of  wave  number   k,   while  the  second  quantity  Is  the 
phase  velocity.   It  states  that  the  speed  of  a  wave  decreases 
with  wave  number  from  the  customary  value  of  unity  (in  our 
normalization) . 

The  same  analysis  Is  now  applied  to  the  full  Navler- 
Stokes  equations,  (1.1-3) •   Taking  transforms  we  obtain 


Initial 


Conditions 


Again  we  are  only  Interested  In  the  dispersion  relation,  now 
given  by 

(1.22)     a^  +  (e  +  6)a^k^  +  (e5k^  +  |  k^)a  +  6k"^  . 

In  keeping  with  the  view  that  the  hydrodynamlcal  equations  are 
valid  for  only  small  wave  numbers,  we  obtain  the  roots  of 
equation  (1.22)  directly  for  small   k. 
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(1.23)  a^~-  I  5k^  +  O(k^) 

(1.24)  a+  ~  t  1/5/3'  k  -  (|  +  |)k2  +  lk^/375^  [  (|  +  |)^+  &£] 

+  0(k^. 


The  mode  o        given  by  (1.23)  Is  strictly  diffusing,  and 
through   6   is  identified  with  heat  conduction.   Equation 

(1.24)  gives  two  propagating  modes,  and  considering  o_^, 
for  instance, 

(1.25)  a^  =  -  (|+  |)k2  +  0(k^ 

(1.26)  F^  =  ^p  -  y^y   (1  +  f  [  (|  +  |)^  +  5£]k^)  +  0(k^  . 

The  first  quantity,   a^,   as  before,  gives  the  attenuation 
rate,  while  (1.26)  gives  the  phase  velocity.   We  now  notice 
however,  that  the  lowest  order  sound  speed  is  the  adiabatic 
value  (in  our  normalization)  and  furthermore  that  the  phase 
velocity  increases  with  wave  number. 

More  accurate  dispersion  relations  may  be  found  by 
choosing  higher  approximations  of  the  Chapman-Enskog  proce- 
dure.  For  instance,  the  Burnett  equations  supply  three 
constants  besides  viscosity  and  heat  conductivity  and  give 
the  roots  correctly  to  O(k^).   Each  higher  approximation 
furnishes  the  next  coefficient  in  a  series  started  by  (1.25-6) 
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Such  investigations  are  to  be  found, for  instance,  in  the  book 

7 
on  ultrasonics  by  Herzfeld  and  Lltovitz.    As  we  have  mentioned, 

we  are  Interested  In  obtaining  a  dispersion  relation  valid 

over  the  full  range  of  wavelengths  and  not  just  an  expansion 

which  furnishes  information  for  large  wavelength  phenomena. 

We  start  this  study  in  the  next  section. 


2 .   Formulation  of  One-Dimensional  Sound  Propagation  According 
to  Kinetic  Models. 

Instead  of  employing  the  Boltzmann  equation  to  describe 
the  molecular  distribution  function,  f,   we  shall  use  a 
kinetic  equation  with  a  mock  collision  term.   The  most 
general  such  equation  to  be  considered  in  this  paper  in, 

dt      ax  V 


mp'/RT^  /RT  ^ 


This  is  the  one-dimensional  (x'),  time  varying  (t')  form  of 

5 
one  of  a  class  of  equations  suggested  by  Gross  and  Jackson"^, 

which  can  be  extracted  from  the  Boltzmann  equation.   For 

reasons  which  will  become  clear,  we  refer  to  this  as  the  triple 
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relaxation  model.   The  following  nomenclature  Is  used: 


^r^  As  /S 


molecular  distribution 


f  =  f(x',tSl^,^2'^3)  '      ^^^ 

function, 

I  =  (I    I  ,  f  )  ,   the  molecular  velocity, 

R,   the  gas  constant, 

p  ,  u  =0,  p  ,    equilibrium  mass  density,  velocity  and 
•^o   o     o 

pressure, 

1 


f     ^ 

r 

P' 

c 

u' 

^ 

f 

1 

p 

J 

_J 

^l/p' 

(C-u-)' 


d^  ,    the  corresponding  non- 
equilibrium  values 


P 


11 


^1  = 


(I  u)  f  d?  -  p',   the  stress  in  the  x-direction, 

— ip (^-u')  f  d^-,,    the  heat  conduction  in  the 

x'  -direction. 


P^ 


■^     2  / 

p'  ^..   e"^^~"'^  ^^■^'^  ,   the  local  Maxwelllan, 


f°  = 


-IV2RT, 


(27rRT  ) 
^     o' 


^n 


,        the  absolute  Maxwelllan, 


m  ,    the  molecular  mass  . 

The  three  symbols   i-l',  V,  y'  >  0   represent  quantities  which 
may  be  obtained  in  the  determination  of  (2.1).   As  given  in 
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reference  5j.  they  are  constants  which  are  obtained  from  the 
molecular  force  law.   In  the  cited  reference  equation  (2.1) 
Is  obtained  for  the  case  of  Maxwell  molecules,   for  which 
the  eigentheory  of  the  Boltzmann  equation  Is  known.    In  fact, 
these  constants  are  merely  related  to  the  eigenvalues  of  the 
linearized  collision  operator. 

For  the  case  of  Maxwell  molecules  equation  (2.1)  has 

sufficient  detail  in  it  to  give  correctly  the  Navier-Stokes , 

4 
Burnett,  and  thirteen  moments  equations.     We  shall  see 

evldeiiCe  of  this  in  later  calculations.   For  other  molecular 

models,  equation  (2.1)  gives  the  values  of  transport  coef- 

o 

flcients  correctly  to  what  Chapman  and  Cowling   call  the 
first  approximation.   Reversing  our  viewpoint,  we  can  pre- 
scribe  [jl',  v',  y'   so  as  to  give  transport  coefficients  or 
combinations  of  them  correctly.   For  Instance,  we  can  pre- 
scribe  i-L '   and   7'   to  furnish  viscosity  and  heat  conductivity 
correctly  (for  a  particular  gas)  and  then  by  a  judicious 
choice  of   v'   approximate  the  Burnett  dispersion  law. 
For  the  case  of  sound  propagation  we  write 


O   .    A 


(2.2)  f  =  f^  +  g 


and  linearize  about  the  equilibrium  given  by   f  .   Linearizing 
the  local  Maxwellian  we  get 


Molecules  which  interact  with  an  inverse  fifth  power  force 
law. 
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(2.3) 


^o-- 


I  Li-  +  p   +  T   V2RT 


O      O       O 


2^     RT 

o 


where 


(2.1+) 


P  = 

A 

T  = 


A 
U 


P' 

T' 
u' 


-  T, 


u 


o 


We  use  the  following  normalization  consistent  with  Section  1 


(2-5) 


■^o 


T 

T  —  m     5 
O 


u 


u 


)^ 


and  in  addition 


( 


P 


'11 


11        p 


S.    = 


^    p^yRT 


o'      o 


(2.6)      <      ^ 


/RT" 


V   = 


£V_ 


m 


o 


^L 


P     U  ' 
'^O 

mv  ' 


T  = 


mv  ' 


(2.7) 


,-^2/2 


cu 


(27r) 


3/2 
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(2.8) 


/  ^  _  1 

T 

V                  1 

T 

)/RT^'        L 

^'.  -    X 

t 

T 

^  =  t 


When  all  this  Is  applied  to  equation  (2.1)  we  get,  finally, 


2 

(2.9)    (|t  +  ^1  ll  +  i)e  =  ^[P  +  (|-  -  |)T  +  ^i" 


+ 


Wii(|  ^1  -  I  ^^)  +  ys  ?.^(|  (^  -   1)] 


The  relation  of  the  "flow  field"  to  the  perturbed  distribu- 
tion function  is  given  by 


(2.10) 


J 


(5  ^^-  1) 


gd^  = 


%    -^ 


u 


^11  ! 

i 


Our  normalizing  constants,   t   and   L,   are  no  longer  the 
mean-free -time  and  the  mean  free  path.   The  procedure 
given  in  reference  5  dictates  the  choice  of   v,   which 
for  equation  (2.1)  is  considerably  larger  than  the  colli- 
sion frequency  of  the  gas. 
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Note  as  In  Section,  we  have  eliminated  the  pressure  with  the 
identity 

(2.11)    p  -  p  +  T  . 

We  again  consider  a  problem  for  which  there  are  no 
boundaries  and  consider  the  initial  value  problem.   In 
Section  6  it  is  shown  that  with  suitable  initial  conditions 
we  can  rigorously  apply  transforms  to  equation  (2.9)-   Operat- 
ing on  equation  (2.9)  with   e         dtdx,   i.e.,  taking  a 
Laplace  transfonn  in  time  and  Fourier  transform  in  space,  we 
get  after  solving  for  the  transformed  perturbed  distribution 
function. 


2 

(2.12)     g  :=  ^ [p  +  (|-  -  |)T  +  i^-n      HP3_^(|  4^  -  i  ^^  ^ 

1  +  a-ik 


+  TS^^^(^  ^^-  1]  + 


5         1+a-ik^^ 

where  we  have  denoted  the  transformed  variable  by  the  same 
letter,  and  given  the  initial  data  a  zero  subscript.   Per- 
forming the  moment  integrations  indicated  in  equation  (2.10), 
and  integrating  out  the  transverse  velocities   ^p,  1^,   we 
get 

(2.13)    V  =  C  V  +  V 


o 


with 
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(2.14) 


V   = 


P 
u 

T 

Pll 
Si 


J 


(2.15) 


V, 


1  +  o-lki. 


(^1    -  ^^/3) 
(^^-  5)^1 


d^ 


(2.16) 


^1 

.2  1 

^1/3     5 

l(^?-l) 


2^(e?-3) 


*^ 


"^1/2 

^2 ^^1        ^■ 

(27r)-^/^(l+a-ik^    )         ^ 


.2      _   1 
'1/2      2 


^^(q-l) 


;2 

'1 


:3 


.2 


7qa  -  I) 


-J         5' 
^i/2-h/2  ^(«i-5i)  i(Si-5«?) 


«?/3-«i/3   ^i/6-^?/3^i    ^{ej-2i?-i)   fe(e^^e^7e,) 


|(?i-5i 


4       _    2fi2      1 
'1/5"    ri~    3 


2^1/6  4 

3 


(q-2^1+2)  ff(^^^^^^l) 


^ 


l(^2_3)l     e5/2-2^^I|,     |(^5_,^3,^^)   l_(^6_,^4^,3^^) 


y 


Each   element    of     C      may  be   evaluated   in   terms    of  a   single 
integral    (see   Appendix   II,    Lemma   l), 
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(2.17) 


M 


A 


T^^^l 


with 


(2.18)     A  = 


1  +  g 
k 


Q 

As  is  well-known  ,  an  integral  of  the  type  given  by  equation 
(2.17)  has  the  real  axis  as  a  natural  cut.   Depending  on 
whether  the  real  part  of  A   is  positive  or  negative  the  inte- 
gral of  (2.17)  defines  two  different  functions  which  are  in 
fact  entire  analytic.    Moreover,  the  analytic  continuation 
of  each  of  these  functions  across  the  cut  does  not  lead  to 
the  other  function. 

At  this  point  of  the  analysis   a     must  be  considered 
as  having  a  positive  real  part,  it  having  bean  so  chosen  to 
define  the  Laplace  transform.   Hence   A  will  have  either  a 
positive  or  negative  real  part  depending  on  the  sign  of  k. 
Equation  (2.17)  thus  defines  two  functions   M   and   M   for 
k  >  0  and   k  <  0,   respectively.   The  precise  forms  of 
these  functions  are  developed  and  given  in  Appendix  I.   It 

is  Interesting  to  consider  the  jump  of  M  across  the  cut. 

9 
This  is  immediately  given  by  a  formula  of  Hermite  ,  and  is 


(2.19)     [M/A]  =  /2T  e   ^' 


a2/2 


where  the  bracket  denotes  the  jump  and  A   the  imaginary  part 
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of  "A.   In  the  Important  limit  when   A.  — >  oo  (corresponding 
to   k  passing  through  zero),  the  jump  is  asymptotically 
small  and  the  transition  between  functions  is  smooth. 

At  this  point, with  all  the  elements  of  the  matrix  computed, 
the  dispersion  relation  is  obtained  by  taking 

(2.20)     det(U-C)  =  0 

where  U  is  the  unit  matrix.  A  study  of  the  complete  dis- 
persion relation  does  not  seem  feasible.  We  will.  Instead, 
make  use  of  expansions,  asymptotic s  and  approximations.  In 
this  way  we  will  map  out  the  dispersion  relation  and  at  worst 
gain  only  qualitative  results.  In  the  next  section  we  will 
be  concerned  with  the  possible  assumptions  that  can  be  made, 
and  their  physical  basis. 


3.   Approximate  Formulations. 

The  work  in  the  dispersion  relation  is  considerably 
reduced  by  the  use  of  the  conservation  equations.   To  obtain 
these,  consider  the  transform  of  (2.9)- 

(3.1)     (a-lk^^  +  l)g  =  co[p  +  (|-  -  f)T  +  i-^u 

+      p(|  ^2  -  1  |2)  ^  ^se_^(|  ^2_  i)j. 
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On  taking  the  first  three  moments  Indicated  by  equation  (2.10) 
we  get  the  transformed  conservation  equations, 


(5.2) 


op    -   iku  =  p. 


(5.3) 


au  -  ik(p  +  T)  -  Ikp-, -,  =  u. 


(3.^4)     aT  -  t  iku  -  I-  IkS  =  T, 
^  3       3        1 


The  subscript  1  refers  to  the  comparable  moment  of  the 
initial  data.  With  the  use  of  equations  (3.2-4),  we  can 
study  the  dispersion  relation  given  by 


(3.5) 


a 

-Ik 

0 

0 

0 

-Ik 

o 

-Ik 

-Ik 

0 

0 

-21k  3 

a 

0 

-2ik/3 

^41 

Ci|2 

C43 

Cj^i^-1 

S5 

^51 

^52 

^53 

C54 

C55-I 

D^(k,a)  =  0 


where  the  C's  are  those  given  in  Table  1.   As  is  clear  by 
comparison  of  (3-5)  and  (2.20),  the  determinant  (3-5)  Is  much 
simpler  to  deal  with.   We  will  therefore  consider  (3.5)  In 
place  of  (2.20).   One  will  notice  that  the  initial  data  prob- 
lem as  dealt  with  in  connection  with  (3-5)  appears  differently 
than  in  connection  with  (2.20).   Later  in  this  section  we  will 
remark  on  the  equivalence  of  the  two  formulations. 
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The  dispersion  relation  as  given  by  (3-5)  will  be  dis- 
cussed in  Section  5-   A  more  complete  discussion  can  be  given 
for  systems  which  are  simpler  than  (3.5).   Instead  of  taking 
(2.16)  as  the  case  to  study  we  can,  for  Instance,  consider 
the  equations  for  (p ,  u,  T)  only,  gotten  by  considering  the 
upper  left  hand  three  by  three  matrix  of  equation  (2.l6),  i.e.. 


(3.6) 


P 


u 


T 


r 


i2v) 


1/2 


;2 

'1/2 


dt 


(1+  o-ik^^) 


J 


X 


.2 
'1 

3 


1/3    3     ^1/3  ^1/3 


p2  1 

^1/2"  2 

^l/2"^l/2 

^1/6-^1/3  ^  ^ 


T 


/" 


+ 


g^d^ 


^ 


1+a-ik^-, 


^y 


(ie"-i): 


.^ 


For  reasons  which  will  become  clear  Ixi  Section  5?  this 
can  be  referred  to  as  the  adiabatlc  case.   This  system  could 
be  obtained  by  setting  U-  =  7  =  0,  or  equivelently ,  by  con- 
sidering instead  of  equation  (2.9),  the  equation 


(3-7) 


4+?iS+l)s  =  -(P+  (^-|)T  +  5,fi) 


or  on  transforming 


(3.8)      (a  -  ik^-L  +  l)g 


^(P  +  (|-  -  f)T  +  e^O)  +  g^ 
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Equation  leads  to  the  dispersion  relation 


(3-9) 


S  = 


^11"^    ^12 


'21 


'31 


'32 


'13 


C22-I   C23 


0,3-1 


=  0 


A  simplification  similar  to  (3-5)  can  be  made  now.   From 
equation  (3-8)  we  see  that  on  taking  the  mass  moment  we 
get  equation  (3-2)  the  continuity  equation.   This,  as  before, 
can  be  used  to  replace  the  mass  moment  In  (3-6).   On  doing 
this  we  get 


(3.10) 


a 

-Ik 

0 

=21 

C22-I 

=23 

=  °3  = 

=   0 

=31 

^32 

S3-I 

which  Is  equivalent  to  considering  (3-9) •   Equation  (3-7)  Is 
referred  to  In  the  literature  as  the  single  relaxation  model 


This  Is  the  name  given  to  equation  (3-7)  by  Gross  and 
Jackson  In  reference  5-   A  relaxation  theory  obviously  has 
all  non-hydrodynamlc  moments  decaying  with  a  single  time 
Inversely  proportional  to   v   of  equation  (2.8). 

We  have  been  somewhat  cavalier  In  setting  [x  and  y   equal 
to  zero,  since  (3-7)  still  Incorrectly  contains   v   In  Its 
normalization.   The  correct  value  as  given  by  Gross  and 
Jackson  Is  related  to  the  collision  frequency  (see  footnote, 
page  18),  which  Is  smaller  than  v   of  (2.8).   We,  therefore, 
think  of  equation  (3-7)  as  having  been  obtained  by  a  process 
In  which  \x     and  y     vanish  and   v   goes  to  the  collision 
frequency.   This  distinction  plays  an  Important  role  In  an 
argument  given  at  the  close  of  Section  4. 
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It  first  appears  in  the  well-known  paper  of  Bhatnager,  Gross 

6 
and  Krook.    In  that  paper  the  authors,  among  other  problems, 

consider  sound  propagation.   Although  they  consider  (3-7) > 

they  restrict  their  calculations  to  a  simpler  model  which  we 

will  consider  in  the  next  paragraph.   Since  our  methods  are 

different  and  since  new  features  are  brought  to  light  we 

will  reconsider  the  case  taken  up  in  reference  4  in  Section  4 

By  further  truncation  of  the  system  (2.l6),  simpler 

systems  can  be  derived.   On  truncating  the  system  to  a  two 

by  two,  we  can  consider 


(3.11) 


'.  "; 


(27r) 


T72 


e"^V2  ^^_ 
1+a-ik^-, 


f. 


^ 


£     ^2 

^1  ^1 


{ 


rj 


/ 


+ 


g^d^ 


1+a-ik^- 


which  is  equivalent  to  choosing  the  model  equation, 

(3.12)  (|^+  ^^  |-+  i)g  =   a3(p  +  ^^u) 
or  the  equivalent  transform 

(3.13)  (a  -  ik^^  +  l)g  =  a)(p  +  i-^u)   +   g^   . 

As  in  the  previous  cases,  one  is  led  to  the  dispersion  rela- 
tion given  by 
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(3.14) 


axi?,- 


l+a-lk^. 


l+a-lki- 


^■^0)  d^^ 


1+a-ikf, 


1 


1+a-ik^ 


=  d; 


0 


Instead  of  the  first  row  of  (3-l4),  once  again  we  can  intro- 
duce the  continuity  equation  which  still  holds,  to  get 


(3.15) 


■ik 


,  ^^jd^^ 


J 


1+a-ik^ 


1 


1+a-ik^. 


D. 


0 


as  the  dispersion  relation.   Using  the  forms  given  in  Table  1 
it  is  straightforward  to  see  that  (3.l4)  is   -  l/l+a   times 
(3.15).   The  previous  two  cases  go  in  the  same  way.   For 
reasons  which  will  become  clear  shortly  this  will  be 
referred  to  as  the  isothermal  case. 

The  most  drastic  truncation  which  can  be  made,  is  to 
consider  the  one  by  one  system,  i.e.. 


.2 
'1/2 


(3.16) 


d^ 


(2Tr) 


T72 


1+a-ik^ 


P  + 


1 


g^d^ 


J 


1+a-ik^. 


which  is  equivalent  to  first  considering 


('•")     4  +  5|j-H)s  =  a^  , 
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or  equlvalently. 


(3.18)     (a  -  Ikl,  +  l)g  =  cup  +  g 


o 


The  dispersion  relation  In  this  Is  simply 


(3.19)    1  -^ 


!   6 


;2 

'1/2 


d? 


2tt 


l+a-lk4- 


1  -  C^^  =  D^  =  0  . 


J 


This  will  be  referred  to  as  the  Isosterlc  case. 

Some  light  can  be  shed  on  the  meaning  of  each  of  the 
cases  In  question  by  studying  the  moment  equations  comparable 
to  the  conservation  equations.   Taking  the  mass,  velocity  and 
temperature  moments  of  (3-7) ?  (3-11)  and  (3-17)?  we  get, 
respectively. 


Sp  ,  Su 


5u  ,  S 


0 


Sp 


(5-20)  \fF^f5?(P+T)  ^W 


11 


0 


^3^3^  ^1~° 


/  ^  +  |H  =  0 
dt    dx 


Bp 


(5-^1)  {  1^  +  It  (p-^)  +  ST 


11 


0 


aT2Su2a      _ 
^^3^"^3^1~ 
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0 


Sp- 


(3-2^)  {  If  ^  It  (p^^)  +  H-" 


=  -u 


So  we  see  that  the  adiabatlc  case  preserves  the  con- 
servation equations,  the  isothermal  case  only  preserves 
the  continuity  and  the  momentum  equations,  and  the  isosteric 
preserves  only  the  continuity  equation  (hence  the  appellation, 
mass-preserving).   The  conservation  equations  (3-20)  should 
be  well  understood.   If  the  stress,  P-i -i  j  and  heat  conduction 
S-,   are  dropped,  the  inviscid  Euler  equations  result.   This 
leads  to  unattenuated  propagation  with  the  adiabatlc  speed 


(=  /5/3   in  our  normalization).   On  the  other  hand  if  the 

Navier-Stokes  relations  for  S   and  p-, -,   are  substituted 

there  results  the  propagation  of  a  diffusing  wave  at  the 

adiabatlc  speed  and  in  addition  a  purely  diffusing  mode. 

A  similar  analysis  can  be  made  for  the  systems  (3-21-23) •   The 

inviscid  theory  is  obtained  by  dropping  quantities  without  a 

time  derivative.   For  (3-21)  this  yields  propagation  with 

speed  one  (the  isothermal  speed,  in  our  normalization),  and 

for  (3-22)  we  get  a  zero  propagation  speed.   In  order  to 

obtain  an  analogous  viscous  description  we  must  develop 

relations  allied  to  the  Navier-Stokes  relations.   To  do  this, 
_  — 

See  reference  10.   Section  6  contains  the  analysis  for  the 
Isothermal  case. 
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we  use  the  method  of  interpolation  which  is  akin  to  the 
Chapman -Enskog  theory,  and  was  invented  by  H.  Grad. 
Essentially,  one  considers  space  derivatives  as  being  smooth 
with  respect  to  time  derivatives.   Under  this  assumption 
the  energy  equation  of  (3-21)  or  (3-22)  can  be  integrated 
(considering  the  space  derivatives  as  constants)  and  after 
sufficient  time  the  solution  is. 


(3.25)     T~-||H.||.s^ 


On  taking  higher  moments  of  equation  (3. 7)  we  can  obtain 
equations  in  p.,-,  and  S-,  (i.e.,  containing  the  time  deriva- 
tive of  these  quantities),  and  again  perform  an  interpolation 
(see  references  h   and  9  foi^  more  details  of  this  method.) 
Neglecting  second  derivatives  we  find  that 


4  Su 
Pll^"  3  ^ 


(3.2i|) 


^12^ 


and  (3-23)  becomes 


(3.25)    T-^-lll 


When  substituted  Into  (3-21)  this  gives 
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t 


+  ^=   0 
dx 


(3.26) 

Su   ,   Sp     o  S  U  _  ^ 


An  asymptotic  study  of  the  fundamental  solution  of  the 
system  (3-24 )  (these  are  the  circumstances  under  which 
(3.24)  can  be  presumed  valid)  shows  that  a  diffusing  wave 
propagates  at  Isothermal  propagation  speed. 

When  the  same  procedure  Is  applied  to  the  Isosterlc 
system  one  gets 


N      p.2 
(3.27)    St  P  -  ^P  ^  ° 


In  this  case,  then,  we  have  diffusion  and  no  propagation. 

The  properties  of  the  various  systems  which  have  been 
sketched  above  will.  In  the  sequel,  be  found  In  a  more 
straightforward  manner.   As  might  be  expected,  the  Isosterlc 
and  Isothermal  cases  allow  the  most  analytical  work  and  are 
the  most  easily  understood.   In  the  following  section  we 
deal  with  these  two  cases  somewhat  exhaustively.   Thereafter, 
the  more  tedious  parts  of  the  analysis  will  be  avoided  and 
only  the  results  stated. 
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4 ,   Isosterlc  and  Isothermal  Propagation. 

In  this  section  we  deal  with  the  two  most  degenerate 
models  discussed  In  the  previous  section.   Both  of  these 
models  are  so  devoid  of  physical  content  that  they  are  useless 
from  this  viewpoint.   However,  many  of  the  properties  of  the 
dispersion  relation  found  In  this  section  will  carry  over 
to  the  more  general  cases  found  In  the  later  sections. 
Because  of  the  nature  of  the  function  M,   an  exact  discus- 
sion of  the  dispersion  law  does  not  seem  possible.   Instead 
we  will  resort  to  a  whole  gamut  of  approximate  devices  which 
will  give  us  a  reasonably  good  description  of  the  dispersion 
law. 

Isosterlc  case. 

Referring  to  (3-19)  we  see  that  the  dispersion  relation 
In  this  case  is  given  by 


(4.1)      k  -  ^  =  F(k,  A)  =  0  . 


As  is  shown  in  Appendix  II,  it  is  sufficient  to  consider 
(4.1)  for  k  >  0  and   Im  A  >  0.   We  will,  therefore,  restrict 
attention  to  the  upper  half  of  the  A-plane  and  use  M   of 
Appendix  I.   Substituting  (l.8)  for  M  into  (4.1)  gives  us 

A 


TT  V2   .2/  1/2 

(^.2)      k  =  (?)     e^  /2  ^^  _  (I) 


J, 
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'  e-^'/2  ^^j 


Also  from  Appendix  I,  we  have  the  asymptotic  expansions 


-  ^  <  arg  A  <  +^ 


and 


(4.4)  M  -,  (2.)V2  e^2/2  ,  1    1   ^  1_2  ,  ...  , 

f  <  arg  A  <  :^   . 

I^tt/^ 
It  Is  Immediate  from  (4.3)  and  (4.4)  that  the  ray  e       ' 

Is  a  Stokes  line  of  M  across  which  the  asymptotic  behavior 

jumps.   This  will  play  an  Important  role  In  what  follows. 

Before  Investigating , the  asymptotlcs  we  examine  (4.2)  for 

real  roots. 

For  real   A   the  right  hand  side  of  (4.2)  is  real. 

From  (4.3)  we  see  that  k  — >  0  as  A  — >•  co,  and  from  (4.4) 

1/2 
k  — >  oo  ad  A  — >  -00.   By  direct  substitution   k  =  (71/2)  ' 

for  A  =  0.   We  see  therefore  that  0  <  k  <  00  goes  into  the 

entire  real  A  line.   If  we  examine  the  roots  in  the  plane 

of  the  Laplace  variable  a,    we  see  that 

for  k  ~  0 

(4.5)  {  a  =  -1     for  k  =  (7r/2)-^/^ 

for  k  — >  00 


-  33  - 


In  what  follows  we  shall  refer  to  the  strip   -1  _  Ra  __  0  , 
and  the  half-plane   RA  >  0  as  the  hydrodynamic  proportions 
of  the  respective  planes. 

We  now  estimate  the  number  of  roots  to  (^.l)  in  the 
A-plane.   To  do  this  we  use  the  method  of  winding  numbers 
(see  [9] J  P-  102).   We  will  consider  closed  regions  only  in 
the  upper  half  plane  by  virtue  of  the  remarks  in  Appendix  II. 
Since  our  expression  (^.l)  is  entire  the  winding  number  will 
only  count  zeros,  and  since  it  is  not  identically  zero  we 

always  choose  a  closed  path  no  point  of  which  is  a  zero. 

12 
Because  of  the  latter  and  Rouche ' s  theorem   we  can  with 

confidence  use  the  asymptotic  evaluations  for  the  determina- 
tion of  the  zeros. 

Let  |A|  be  large  and  consider  the  closed  path  of  the 
diagram.  Figure  4.1.   For   A   sufficiently  large,  the  leading 
terms  of  the  dispersion  relation  are 


(4.6)      0  =  k  -  (27r)^/^e'^  /^  -  ^  +  0  (l/A^) 


2  <  arg  A  <  TT 


(4.7)      0  =  k-Y+  0(1/^^)  0  <  arg  A  ,  7r/2 


A 


The  angular  splitting  which  is  somewhat  arbitrary  is  made 
merely  for  convenience  in  reference.   We  see  that  neither  (4.6) 
nor  (4.7)  (or  for  that  matter  (4.1))  can  have  roots  in  R  as 
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k  — >  oo. 

Next  consider  the  region  of  Figure  2.   To  consider  this 
region  we  make  use  of  the  fact  that  if  Y^.  ■]  F(A,k)j  vanishes 

at  2n  distinct  points  on  a  +  b  +  c ,   F  has  at  most  n  zeros 

12  I  I 

in  R.    For  any  fixed  value  of  k  as  1^1  — >  oo,   we  see 

directly  from  (4.2)  that  1^  [F]   has  at  most  one  zero  on 

a  +  b,   and  from  (4.3)  that  it  has  no  zero  on   c.   Hence, 

with  the  possible  exception  of  the  zero  on  real  line, 

already  discussed,   F(k,A)  can  have  no  zero  in  the  right 

half  of  the  A-plane. 

~^TT  TV 

We  now  consider  the  region  of  -4^  >  arg  -^  >  o"-  P^om 
the  asymptotic  expansion  given  by  (4.4),  we  see  that  for 
any  |k|  >  0,  no  matter  how  small,   (4.2)  will  not  have  a 
solution  for  sufficiently  large  distances  from  the  origin 
along  a  ray.   This  is  true  for  any  ray  in  the  region  as 
close  as  we  please  to  the  Stokes  line.   Of  course,  the  closer 
such  a  ray  is  to  the  stokes  line,  the  further  out  we  must  go 
to  guarantee  non-solvability.   One  must  not  conclude  that 
the  solutions  are  on  the  Stokes  line,  but  rather  that  they 
are  closer  than  any  ray  (for  instance  a  parabola  is  closer 
than  a  ray).   We  now  further  explore  this  region,  i.e., 
region  I  given  in  Figure  3.   We  again  make  use  of  the  theorm 
on  the  ^{F(k,A)j.   From  (4.1)  and  (4.6)  for  a  fixed  k  and 
|A|  sufficiently  large,   I^  F  has  no  zero  on  b,   and  at 
most  one  zero  on  c.   Along  the  Stokes  one  may  easily 
estimate  the  magnitude  of  the  right  hand  side  of  (4.2)  to 

4 
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be  <  (7r/2)    +  /2"   (asymptotically  the  estimate  is  /27r 

which  is  slightly  smaller).   Therefore,  the   I^'iFJ  will  not 

1/2    — 
have  zeros  along  the  Stokes  line  if   k  >  (7r/2)  ^   +  /2  . 

And  with  the  proceeding  information,   F  will  not  have  zeros 

in  I  for   k  >  {ii/2)^'      +  /2    .      We  will  later  asymptotically 

locate  these  roots. 

To  estimate  the  number  of  zeros  in  region  II  of  Figure 

3,  we  compute  the  winding  number  in  traversing  the  path 

a  +  d  +  e.   The  only  contribution  will  come  from  a  +  d. 

First  we  consider  the  arc  d.   The  winding  number  is  obtained 

from 


(4.8) 


1,-1 
TT-  tan 
27r 


;27r) 


1/2  \i^-v^/2       . 
'       e  '       sm  i-LV 


k-(27r)l/2  ^\?-v^/2    ^^3  ^^ 


where 


(^•9) 


A  =  |j.  +  iv 


=  Re" 


on  d . 


On  the  arc  a  the  product  |i.v  decreases  from  R  /2  to  0. 
Using  only  asymptotic s  on  the  path  a,  the  winding  number 
is  gotten  from 


(4.10) 


1 
27r 


tan 


-1 


(27r)^^^sin  p^ 
k-(27r)^/^  cos  p^ 


where 
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(J+.ll)    A  =  p(-l  +  i)  . 

The  Important  thing  to  note  Is  that  if 

(4.12)    k  <  (27r)^/^ 

the  contribution  from  (4.10)  "unwinds"  the  contribution 
from  (4.8).  In  this  case  the  number  of  roots  Is  0(l). 
If, on  the  other  hand. 


(4.15)     k  >  (27r)l/2  ^ 


(4.10)  does  not  contribute  and  It  Is  easily  seen  from  (4.8) 
that  the  number  of  roots  In  II  Is  0(R  ).   One  notices  In 

(4.8)  that  the  denominator  will  not,  for  a  fixed  large  k, 

2    2         2 
vanish  unless  m-  -  v   =  0(^nk  ).   This  Is  In  agreement  with 

our  previous  remark  on  the  decrease  of  the  number  of  zeros 

In  a  bounded  region  as  k  grows.   We  can.  In  fact,  very 

easily  calculate  from  (4.8)  the  "exact"  number  of  zeros  In 

II  for  k  >  (27r)"'"/  .   We  will  not  do  this  but  rather  go  on  to 

the  location  of  the  zeros. 

We  first  asymptotically  locate  these  roots  which  lie  on 

the  Stokes  line.   On  the  Stokes  line  we  may  write,  on  using 

(4.11)  In  (4.6),  and  splitting  Into  real  and  Imaginary  parts 

(4.14)     k  =  (27r)^/^  cos  p^  +  ^ 
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(4.15)     (27r)^/^siri  p^  = 


^)^ 


For 


(4.16) 


I  p  ~  (27rN)  -  e 


A  I  '-2(7rN)^/2 


equation  (4.15)  gives 


(4.17)    £ 


2Tr/N2 


and  from  (4.l4) 


k  --  (27r) 


1/2 


For  further  reference  It  Is  worth  noting  that 


^  p^^  (2N+l)7r  +     ^ 


(4.18) 


27r/2N+l 


A  I  ~  /2(2N+l)7r 


1/2 


k  ^   -(27r) 


1/2 


On  returning  to  (4.2)  we  get  on  differentiation  that 


(4.19) 


dA  ^  _1 

dk   Ak-1 


Using  the  notation  of  (4.9)  we  find  the  following  equations 
for  the  real  and  Imaginary  parts  of  (4.19). 


(4.20) 


d[i- 
dk 


^ik  -  1 


\^{\i^+   v^)-2^ik  +  1 
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(4.21) 


dv  ^  ^vk 

'^^        k^i[i^+   v^)-2\ik   +  1 


Taking  the  ratio  of  the  two  equations, 


^^"^"^1        d^i   lik-i 


From  our  previous  estimates  we  know  that  in  the  neighborhood 
of  the  Stokes  line  and  to  the  left  of  it   M-k  »  1.   Under 
this  assumption. 


(4.23) 


dv  ^  -V 
dix  ~  \i 


and  on  integrating 


(4.24)     V 


-7rN2 
1^ 


The  information  along  the  Stokes  line  given  by  (4,l6)  and 
(4.18)  has  been  used  as  initial  data  to  get  the  solution 
(4.24).   This  information  along  the  Stokes  line  also  states, 
by  virtue  of  (4.23),  that  the  branches  of  the  dispersion 
relation  are  all  normal  to  the  Stokes  line.   Further,  from 
(4.20)  and  (4.21)  we  have  that 


(   ^ 
dk 


^  =  0 


(^.25)  < 


dM- 
dk 


^V^     _      _-i 


for  k 


0 
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and 


(4.26) 


dk 


dv 
dk 


=  0 


-1 
vk 


for  k  =  -, 


note   k  <  0 


With  the  information  gathered  in  the  preceeding,  we  may, 
with  reasonable  accuracy,  sketch  the  dispersion  curves  in 
the  A-plane.   A  typical  curve  is  plotted  in  Figure  4. 
Actually,  of  course,  an  infinite  number  of  curves  like  the 
one  shown  should  be  plotted.   The  dotted  curve  in  the  dia- 
gram indicates  the  value  for  k  <  0  which  does  not  belong 
to  the  dispersion  relation,  and  is  plotted  only  for 
completeness.   One  should  also  notice  that  the  real  axis 
is  also  part  of  the  dispersion  relation.   In  fact,  the 
contribution  in  the  hydrodynamic  half-plane  is  asymptotically 
the  most  important  part  of  the  spectrum.   In  connection  with 
this  plot  we  make  one  last  calculation.   For  the  region  to 
the  left  of  the  Stokes  line  we  can  neglect   v   compared  to  [i 
and  if  (4.20)  is  solved  under  this  assumption,  we  get 


(4.25) 


^^ 


in  ^  +  27rN 


and  from  (4.24) 


(4.26) 


ttN 


£n^  +   2-vrN 
27r 
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since 


(^.27) 


a  =  kA  -  1 


we  can  give  the  sketch  (Figure  5)  of  the  spectrum  In  the 
a-plane . 

Isothermal  case 


It  is  not  our  Intention  to  go  through  the  same  for  this 
case  as  In  the  Isosterlc  case  just  considered.   With  the 
exception  of  the  hydrodynamlc  part  of  the  plane,  the  same 
features  of  the  dispersion  equation  are  found.   For  this 
reason,  only  a  sketchy  analysis  will  be  given.   The  hydro- 
dynamic  part  of  the  plane, on  the  other  hand,  will  receive 
more  care. 

The  determinant  in  this  case,   Dp,   Is  given  by  (5-15), 
and  in  a  slightly  different  form  is 


(^.27) 


D. 


kA 


I  (1-M) 


-Ik 


£  (l-M)-l 


=  0 


On  expansion,  the  dispersion  relation  is  given  by 


(4.28) 


0   =  -D^   =.    (M-1)  +  A^(l-M)  +  1  -  ^  (1-M)  -  kA 


We  first  restrict  attention  to  the  second  quadrant  of  the 
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A-plane.   In  this  region  the  asymptotic  expansion  is  given 
in  Appendix  I   as 


(4.29)         M  -  (27r)^/2  Ae^  /^  +  1  -  ^  +  ^  -  . 


With  the  experience  gained  in  the  isosteric  case,  a  few 
short  calculations  give  us  the  dispersion  law  in  this 
region.   First,  we  consider  the  asymptotic  roots  for  k  =  0, 
The  dispersion  relation  then  is 


P 
(4.30)  (27r)l/2  ^A  /2^  1^   Q(l  ) 


This  implies  that  these  roots  all  lie  to  the  right  of  the 

Stokes  line.   To  see  this  we  recall  that  along  the  Stokes 

A^/2 
line   e  '^   is   0(l)   whereas  along  any  ray,  to  the  right 

is  exponentially  decreasing  and  any  ray  to  the  left,  exponen- 
tially increasing.   Next  it  is  useful  to  find  the  roots 
along  the  Stokes  line  as  A  — >  co.   it  is  simple  to  see  that 
the  only  possible  k's  for  which  there  are  roots  are  those 
for  which   k  —7-  00.   in  this  case,  then,  the  dispersion 
relation  is  approximately  given  by 


(4.31)         -D    ^   A(l-M)  -  k  ^  0 


Introducing  (4.11)  we  may  write  the  dispersion  relation  as. 
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(^.32)      ,   21p2  e-^P^  =    \.^   +   0(1) 


(27r)- 


This  leads  to 


fh    ^x^  ^2    (2N-1) 

(^•33)        p  • —  ^^ — Q— ^  IT 


2 


and 


(^.3^)         ^ttp  -^  t  (2N-l)7r 


Finally,  we  consider  the  asymptotics  in  the  region 
^  <  arg  A  <  TT.   Looking  at  (^.28)  the  same  terms  are 
dominant ,  i.e.. 


(4.35)       ■  D5^A^(i-M)  -  kA— -0 


This  being  the  case,  each  branch  may  be  associated  with 
the  point  at  which  it  crosses  the  Stokes  line.   We  will  be 
content  with  just  the  asymptotic  estimate  of  A  which  a 
simple  calculation  shows  to  be 


(4.36)         A-Vin  1^  -  (N-l)i7r 


The  three  short  calculations  just  made  tell  us  that  the  left 
half  of  the  A-plane  would  have  a  form  comparable  to  that 
given  in  Diagram  4.   We  will  sketch  this  shortly,  after 
having  considered  the  hydrodynamic  region. 
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Rather  than  perform  the  appropriate  asymptotlcs  with 
the  expansion  of  Dp,   it  is  more  revealing  to  perform  the 
asymptotlcs  directly  in  (4.27).   In  this  form  the  asymptotlcs 
give 


(4.37)    D. 


-Ik 


Ik 


(1+a) 


^ 


i+g  / 


\e 


3k 


k 


(l+a)2    (i+a) 


T 


•)  -  1 


0 


For  convenience  we  have  returned  to  the  Laplace  variable  a. 
Performing  row  operations  this  can  be  put  into  the  form 


(4.38) 


Ik 
1+a 


-Ik 


3k' 


(1+a)' 


=  0 


Inspection  shows  that  the  only  roots  of  (4.38)  consistent  with 
A  — ^  00  are  of  the  form  a  =  0(k).   Imposing  this  condition 
and  again  making  use  of  row  operation,  we  get 


(^.39) 


-Ik 


-Ik 


a+  2k 


=  0 


On  comparing  this  with  the  system  (3-26)  we  see  that  it  is 
precisely  the  Fourier-Laplace  transform  of  that  system.   The 
Chapman -Enskog  procedure  performed  on  that  system  is  equivalent 
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2 
to  an  expansion  In  small  wave  number.   To  0(k  ),  (4.39)  gives 


(J|.40)         a  =  t  Ik  -  k^ 


Therefore,  to  this  order  we  get  two  attenuating  waves, 
traveling  with  the  isothermal  speed,  to  the  right  and  left 
Translating  (^•.■H-O)  Into  the  A-plane  we  have 


(4.41)     11m  A  -^  +»  t  1 
k— >0 


Collecting  together  the  calculations  made  for  the  isothermal 
case  we  have  Figure  6,  the  A-plane. 

Representing  the  dispersion  relation  in  the  a-plane  we 
get  Figure  7-   On  comparison  with  Figure  5  we  see  that  the 
primary  difference  is  the  propagating  modes  which  occur  in 

the  hydrodynamic  portions  of  the  plane,  for  the  isothermal 

13 
case.   Using  existing  tables   ,  a  numerical  calculation 

shows  that  the  hydrodynamic al  branches  leave  the  hydro- 
dynamical  portion  of  the  plane  at 


(4.42) 


T'        ,  when  the  absorption  time  is  equal  to  the  collision 
time  ti   .  I  is  approximately  lOfo   over  its  isothermal 
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value . 

In  both  cases  studied  In  this  section  we  encountered 
a  peculiar  behavior  of  the  dispersion  relation  In  the  region 
I^  a  <  -1.   Since  the  frequencies  In  this  region  exceed  the 
collision  frequency,  we  shall  refer  to  this  as  the  Knudsen 
spectrum.   It  is  indeed  very  closely  related  to  free-flow. 
Referring  to  equation  (3-7)  for  Instance,  we  see  that  if 
the  dispersion  law  is  plotted  in  a  dimensional  a-plane  and 
V  — >  0,   that  the  hydromagnetlc  region  narrows  to  a  thinner 
and  thinner  slit.   As  we  approach  the  limit,  the  "Knudsen" 
spectrum  becomes  the  left  half -plane.   The  limit   v  =  0, 
however,  cannot  be  assumed  since  the  dispersion  relation  is 
wiped  out  identically.   This  state  of  affairs  is  directly 
due  to  a  certain  property  of  the  model  equations  we  are 
studying.   It  is  shown  in  reference  4  (see  also  Section  6) 
that  solutions  of  such  equations  evolve  in  time  according  to 
a  finite  number  of  moments.   For  Instance,  the  solution  to 
equation  (3-7)  Is  time-dependent  only  through  its  functional 
dependence  on  p ,  T,  and  u.   For  free  flow  there  is  no  basis 
for  such  a  functional  dependence.   This  accounts  for  the  singu- 
lar behavior  of  the  dispersion  relation  when   v  =  0. 

On  referring  back  to  equation  (2.17)  defining  M,   we 
see  that  the  Knudsen  spectrum  occurs  in  connection  with  the 
analytic  continuation  of  the  basic  Integral,   M.   This  same 
effect,  i.e.,  of  a  fountain  of  lines  from  a  =  -1,   occurs 
for  the  single  relaxation  model  (3-7)?  for  our  canonical 

-  46  - 


five  moraents  case  (2.1),  and  Indeed  for  all  generalizations 
of  (2.1)  given  by  Gross  and  Jackson.   The  "fountain"  occurs 
specifically  at  a  =  -1   only  because  of  our  normalization, 
and  In  a  dimensional  framework  would  occur  at  a  =  -v . 
Recalling  the  remarks  made  In  connection  with  footnote  9, 
we  see  that  the  location  of  the  fountain  of  times  for  the 
single  relaxation  model  occurs  to  the  right  of  a  =  -v, 
say  at  a  =  -v  .   In  the  silt  region  of  the  [x-plane  given  by 
-V  <  R  a  <  -V   the  five  moment  system  gives  something 
entirely  different  than  the  fountain  predicted  by  the  single 
relaxation  model  In  this  case  (see  next  section).   A  better 
approximation  than  the  five  moments  case  would,  In  turn, 
move  the  location  of  the  fountain  further  to  the  left  In 
the  0-plarie.   It  Is,  In  fact,  the  case  that  the  exact  disper- 
sion relation  does  not  contain  the  fountain  nor  anything 
resembling  it. 

The  fountain  1b  only  a  peculiarity  oi'  our  equation  and 
has  no  basis  In  physical  reality.   l<'or  this  reason  It  will 
not  be  further  considered.   On  thi   iin'r  hand,  the  hydro- 
dynamic  specti'um  has  many  reassuring  features. 
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5-   Single  and  Triple  Relaxation  Equations. 

Single  Relaxation  Dispersion  Law 

Continuing  the  program  outlined  in  Section  3?  we  now 
look  into  the  single  relaxation  model  given  by  equation  (3-7) 
The  dispersion  relation  for  this  case  is  given  by  (5-10), 
which  upon  simple  row  and  column  operations  becomes. 


(5.1) 


^3  = 


kA  - 


1  ■ 

k 

"   A 

-ik 

0 

1 

A 

C22-I 

C23 

0 

^32 

C33-I 

=  0 


On  expanding  the  determinant  we  may  write 


(5.2) 


k^  -  Z^k^  +  Zgk  -  Z^  =  0 


where 


(5.3) 


„     A^(M-l)  ,  IIM    2AM  ^  5^ 
\   =     6    +  SA ^  +  ^ 


(5.4) 


_  2M^    2M(M-1)  ,  A^(M-l)  _  2M 
^2-^-     3     ^     3     "3"''' 


(5.5) 


„   _  -2M(M-1)  ,  A(M-l)  ,   M 
^3  ~    3A    ^  ~b  ^  EX 


We  are  first  interested  in  the  asymptotic  roots  corresponding 


to  k  '^  0.   This  will  correspond  to  | arg  A|  < 


37r 


Rather 
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than  performing  the  asymptotlcs  directly  In  equation  (5-2) 
it  Is  more  revealing  to  perform  the  asymptotlcs  directly  In 
the  determinant,  (5.1).   On  returning  to  the  {o ,k)    notation 
this  Is, 


(5.6) 


■Ik 


0 


-Ik 


4  ,2 
a+  ^  k 


5 


Ik 


■Ik 


0 


In  order  to  obtain  this  we  have  made  use  of  the  fact  that 
all  roots  are  at  least  such  that   a  =  0(k).   Furthermore, 
the  entries  have  been  calculated  to  O(k^).   Comparison  of 
(5.6)  with  eauatlons  (3.20)  and  (3.24)  shows  that  (5.6) 
can  be  Identified  with  the  Pourler-Laplace  transform  of 
that  system.   This,  as  mentioned  before,  Is  formally  the 
linearized  Navler-Stokes  equations.   Equations  (3-24)  show 
that  the  ratio  of  viscosity  to  heat  conductivity  Is  off  a 
factor  of  2/3  from  the  first  approximation  of  this  ratio  as 

Q 

given  by  Chapman  and  Cowling  .   This  has  already  been  noted 
In  the  literature."*"  '"""^  That  this  should  be  so  Is  evident 
from  the  fact  that  the  single  relaxation  model  contains 
only  a  single  constant  (which  we  used  to  normalize  the 
space  and  time  coordinates).   As  pointed  out  previously, 
this  single  constant  can  be  adjusted  to  give  either  the  vis- 
cosity or  heat  conductivity  correctly,  but  not  both,  in 
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general.   The  presence  of  only  one  constant  prevents  the 
results  of  the  single  relaxation  model  from  being  anything 
but  qualitatively  correct.   For,  in  order  to  have  agreement 
for  small  wave  number  phenomena  with  the  Navier-Stokes 
theory,  two  constants  are  necessary.   The  three  relaxation 
model,  taken  up  later  in  this  section,  does  fulfill  this 
minimum  requirement. 

The  determination  of  the  roots  of  (5-l)  in  the  asymp- 
totic limit  is  gotten  in  a  straightforward  manner  and  one 
obtains. 


(5.7)  a  =  -k^  +  O(k^) 


2 
(5.8)  a   =  ±   '/VI  ik  [1  +  ^]  -  k^  +  O(k^) 


The  first  root  (5-7)  is  pure  diffusion  and  is  associated  with 
hydrodynamic  heat  conduction.   The  two  roots  given  by  (5-8) 
involve  propagation  to  the  right  and  left  as  well  as  diffti- 
sion.   One  notices  that  to  0(k)  the  propagation  speed  is 
given  by  Y'^/ 3  ,  which  is  the  adiabatic  speed  in  our 
normalisation.   (For  this  reason  we  sometimes  refer  to  this 
as  the  adiabatic  model.)   To  0(k  )  diffusion  enters  and  to 
0(k  )  we  obtain  an  increase  in  the  phase  speed. 

Using  the  tables  of  reference  15,  a  calculation  of  the 
phase  speed  when   I^  a  =  -1   was  made.   The  results  are 
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Im  ^  -^  1.307 


(5.9) 


k 


1.35 


One  sees  that  the  phase  speed  is  only  slightly  larger  than 
the  adlabatlc  speed  ( — •  1.307)5  of  the  order  of  one  per 
cent.   In  view  of  the  relation  for  k  '^^  0   given  by  (5-8), 
this  indicates  the  interesting  possibility  of  the  phase 
speed  achieving  a  finite  maximum  for  0  <  k  <  1.35-   In  a 
report  to  follow,  a  full  numerical  analysis  of  this  situa- 
tion will  be  given. 

Using  the  above  calculations  we  get  Figure  5-1-   A 
sketch  of  the  "fountain"  has  been  included  since,  as  inspec- 
tion shows,  it  is  slightly  different  than  the  isothermal 
case  given  in  Figure  4.7- 

Triple  Relaxation  Dispersion  Law. 

As  shown  in  Section  3>  rather  than  consider  equation 
(2.20)  we  can  study  the  three  relaxation  time  system  with  the 
conservation  equations  already  included. 


(5.10)    D^  - 


a 

-ik 

0 

0 

0 

-ik 

a 

-ik 

-ik 

0 

0 

-2/3   ik 

0 

0 

-2/3   Ik 

^41 
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C44-I 

C45 
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^52 

=  55 

C54 

=55-1 
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As  usual  the  last  two  rows  are  to  be  evaluated  from  Table  I 
If  the  determinant  is  expanded  and  the  entries  substituted 
from  Table  I,  we  get  the  following  dispersion  relation: 


(5.11) 


where 


(5.12: 


0 


5   ^ 


k5-  k^Z^  +  I^+  E^L]  +  k5[Z2-  n[|(l-  i 


2 


+  2B(M+  1  +  1^)]  +||a  +  D(^-  5)j 


+.T  |(A2+f)  ]  -k2[Z3+.   §^.^ 


+  J|^+|g^,+ 


10 


15^ 


15 


9A^-  5|]^,[HiMB^g,M^] 


5A 


10k  C 


A  = 


B 


55 


yA 


^  ^^^41 


_A[  -2?\^  +  5A^r 


-  M^(5  +  3A^) 


V 


D 


E  = 


hV?    -    (A^+  1)M  +  A^  +  6 


The  calculation  of  (5. 11)  is  tedious  and  of  a  type  that  one 
does  not  like  to  often  repeat.   It  is,  furthermore,  virtually 
impossible  to  analyse  in  its  entirety.   However,  a  report  to 
follow  will  contain  a  numerical  analysis  of  (5-11) • 
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The  remainder  of  this  section  will  be  devoted  to  the 
asymptotic  study  of  the  triple  relaxation  time  model.   As 
one  can  see,  this  relation  is  very  cumbersome  and  difficult 
to  analyse.   We  will  only  consider  the  asymptotics  in  this 
section.   We  will  also  make  plausible  conjectures  on  the 
remainder  of  the  spectrum  based  on  what  we  have  found  in 
the  previous  sections.   In  a  report  to  follow  we  will  give 
a  full  discussion  of  D^   from  numerical  data.   As  before, 
rather  than  considering  the  expansion  of  D^-,   we  will 
perform  the  asymptotics  directly  in  (5-10).   Using  the  forms 
of  Table  I,  we  get 


(5.15)   0  -  D, 
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The  only  roots  of  (6.2)  consistent  with  the  asymptotic 
expansion  are  those  for  which  k  — >  0.   Further,   0(l/A  )  is 
only  retained  in  the  first  three  columns  since  expansion  shows 
that  those  appearing  in  the  fourth  and  fifth  columns  contribute 
orders  which  have  already  been  neglected.   Another  equivalent 
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way  of  seeing  this  Is  in  noticing  that  if  (5- 13)  is  considered 
as  the  matrix  of  the  system,  then  the  fourth  and  fifth  columns 
multiply  the  stress  and  heat  conduction,  respectively.   In  the 
limit  of  small   k,   both  these  quantities  are  0(k)  whereas 
the  density,  temperature  and  velocity  are  0(l).   Replacing  A 
in  favor  of  a  and  k,  we  reduce  (5-13)  to, 


(5-1^)   0  =  Dj 
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■ik 
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Reduction  to  (5.l4)  is  tedious  to  explain.   In  short,  it  is 
gotten  by  multiplying  a  row  by  the  proper  coefficient  and 
adding  it  to  another  row,  always  neglecting  higher  orders. 
Setting   k  =  0  in  (5.l4),  we  get  the  roots 


(5.15) 


0   =   0,    -1+2^1,  -1+Y 


where  the  zero  is  counted  thrice.  The  triple  root  corres- 
ponds to  the  three  branches  found  in  the  single  relaxation 
model.   The  other  two  branches  implied  by  the  non-zero  roots 
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will,  as  we  shall  see  shortly,  Introduce  some  novel  effects. 
We  first  deal  with  the  triple  branch  emanating  from  the  origin, 
and  which  we  will  call  the  hydrodynamic  branch. 

In  this  case  as  in  the  previous  treatment  of  the  hydro- 
dynamic  part  of  the  plane  we  may  write 


(5.16) 


2      ^ 

a  =  ak  +  pk  +  yk  + 


and  then  evaluate  the  constants  a,  p,  y,  ...  .   Before  doing 
this  we  observe  that  for  these  roots,  since  a  =  0(k),  we  may 
put  (5.l4)  in  the  following  form 


(5.17)   0  =  D, 
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This  is  precisely  the  form  of  the  dispersion  relation  which  is 
gotten  from  the  thirteen  moment  equations  in  our  normalization, 
(see  reference  4).   However,  it  can  only  be  applied  to  the 
three  branches  which  have  the  form  given  by  (5.16).   As  we 
will  see  shortly  the  remaining  roots  behave   quite  differ- 
ently than  what  is  predicted  by  the  thirteen  moments  equations. 
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Returning  to  equation  (5.l6)  we  have,  on  substituting  into 
(5-17)^  the  three  branches 


(5.18)  o  =  |rT+  o(k'') 

(5.19)  .-txrmi^^v^   (5(r-l)(2.-i)  ^  7^^ 

'  15(2'-!)^"  "  •^"'irFTT  ^  K2f::TT  '  ^  oC^*' 


Comparison  with  the  similar  results  found  for  the  single 
relaxation  model  shows  that  the  expansions  disagree  at  the 
0(k  )  term  and  onward.   Of  course,  as  was  pointed  out  in  the 
discussion  of  the  single  relaxation  model,  only  one  constant 
is  given  in  that  theory.   Thus,  for  instance,  we  can  only 
give  either  viscosity  or  heat  conductivity  correctly.   On 
the  other  hand,  since  the  three  relaxation  time  theory  has 
three  constants  built  in,  it  is  capable  of  more  latitude. 
A  study  of  the  exact  linearized  Boltzmann  equation  made  in 
reference  4  shows  that  (5-l8)  and  (5-19)  give  the  correct 
result  to  0(k-^)  . 

In  examining  the  other  two  branches  indicated  by  equation 
(5.15)  we  can  no  longer  avail  ourselves  of  (5.17)  since  the 
expansion  used  in  deriving  it  is  no  longer  valid.   This  then 
tells  us  that  the  behavior  in  the  neighborhood  of  the  two 
remaining  branches  will  differ  from  that  predicted  by  the 
thirteen  moments  dispersion  relation.   A  straightforward 

-  56  - 


calculation  gives  for  these  two  branches 


2 


(5.20)    a  =  -l-r+k^[l5^.|(^2f;3fr^)]+0(k*) 


If  the  values  for  a  Maxwell  gas  for  \i     and  y     are  used, 
one  gets  the  sketch  given  In  Figure  5-2.   By  slightly  varying 
the  parameters  \i     and  y ,      we  can  get  the  picture  sketched 
In  Figure  5- 3. 

A  detailed  analysis  of  the  exact  Boltzmann  equation, 
however,  shows  that  the  circled  regions  of  both  diagrams  of 
this  section  are  incorrect.   For  a  detailed  discussion  of 
this  point  the  reader  is  referred  to  [4].   The  latter  report 
also  shows  that  the  hydrodynamic  branch  is  given  exactly  to 
0(k  )  in  the  neighborhood  of  origin.   We  further  expect  it 
to  give  at  least  a  good  qualitative  description  in  the  remainder 
of  the  region. 
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6.   General  Properties  of  Model  Equations. 

The  model  equations  which  we  have  been  considering  have 
been  systematically  developed  from  the  Boltzmann  equation  by 
Gross  and  Jackson  [5]-   Although  their  procedure  is  quite 
elegant,  it  is  without  rigorous  foundation  and  from  the 
mathematical  viewpoint  is  strictly  ad  hoc.   Further,  the 
physical  basis  for  any  model  equation  is  at  best  obscure. 
This  being  the  case,  one  would  like  some  assurance  of  the 
reasonableness  of  problems  formulated  with  the  model  equa- 
tions.  This  section,  which  stands  apart  from  the  rest  of  the 
report  in  spirit  (although  it  depends  upon  the  previous 
results),  will  demonstrate  some  essential  features  of  the 
initial  value  problem  for  the  model  equations,  namely  that 
solutions  exist  and  are  unique  and  that  they  asymptotically 
approach  hydrodynamics. 

Existence  and  Uniqueness  of  Solutions 

We  investigate  the  following  equation 

A 
>_ 

n=0 


(6.1)  l#  +  ^  1^  +  g  -  ^  X~  aA^ 

^         ^  6t  ax        ^  ^ — 7T   n  n^n 


where 


,'   5mn   J 


(6.2) 


00  i)   d^ 

—  00 


an  =  j  g^^d-^  : 


—  oo 
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The  X  are  constants  and  the  ^  =  Ai^^^)  polynomials  in  ^ 
For  the  following  considerations  the  sign  of  A  Is  Im- 
material.  As  before  we  denote  the  Gaussian  by 

(uo)  60  = 


(27r) 


ys 


One  should  further  note  that  we  are  no  longer  restricting 
attention  to  one  dimensional  phenomena.   It  Is  clear  that 
all  models  mentioned  earlier  are  of  the  form  of  equation 
(p.l).   We  denote  the  initial  conditions  of  (6.1)  by 


(6.4)  g(x,^,0)  =  g^ 


Rather  than  consider  the  quality  g,   we  equivalently 
consider 

(6.5)         G  =  4^- 

Using  an  obvious  notation  we  have 

N 

(o.b)  ^+  j.  ^  =  y  f,  ^  ^ 

^         '  at  dx       "'^     n   n-^n 

n=l 


(o.V)  G(t=0)    =  G 
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(6.8)  A^   =   /  coGjd^di 
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Regarding  the-   A,  {=   A    (/, ,t))      arj   known   functlonrj,    v/o   ma.}/ 
Integrate    (6.6)    to   obtain, 


r*    " 


(6-9) 


with 


G  =  Gjy.-<=.t,    i)   + 


'>         A    {y.-'^/-,,v/)r  /j  Chi 
"^—^r       n  n^  ri 


o 


n=l 


(6.10) 


X     =  z-i^[t-r3] 


It  Is  easily  seen  that  solving  (6.^')  Is  equivalent  to 
solving  (6.6).   We  shall  shov/  existence  and  ijnlqueness  of 
solutions  to  (6.9)  In  the  usual  v/ay.  I.e.,  by  shov/lng  that 
the  following  Iteration 


/  II 


(6.11) 


G 


^''■"^^  =  G,_^(z-tt,e)  +    X  4''^(^*)Vk''^ 


k=l 


defines  a  contraction  mapping.  To  shov/  this  v/e  v/111  assume 
cert8-ln  conditions  on  a  finite  number  of  moments  of  the 
initial  d3.ta.   V/e  assume  the  existence  of 


(6.12) 


\iO  "  ^  ''^^-  ''^-''  ^r,(^-'^0''^^^ 


•n  -3;-  o 


Actuall?/  instead  of  (6.11)  v/e  consider  equivalently. 


# 


It  is,  in  fact,  shovm  in  the  cited  reference  tl-i^-t  any 
finer  model  v?ill  not  alter  the  loca,tion  of  the  roots  to 
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(6.13)    A^'^+l)  = 
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\ 
t  ,   N 


-i 


j|^A(")(x")A^^^^a3dl  ds 


This  should  be  regarded  as  a  matrix  equation 


(6.14) 


a("+i) 


K- 
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M  ,  Af^'^ds  +  A  , 
M-k   k        oi-L 


where 


(6.15)    M  ,  =  [/  <A^   j^,  CD  A,  ] 
^     '  |J,k    •'    '^M.'^k   k-' 


and  Is  an  operator,  and 


(6.16)    \  =  [/  ^Gj^^^d^]  , 


k  =  1,  .  .  .  ,  N 


and  similarly  for  A  ,  .   As  a  norm  we  take 


(6.17)    ||A(x,t) 


max[  A-,  (x,t )  ,  A5(x,t), 
x,t   ^        ^ 


.  A^(x,t)], 


I.e.,  we  take  the  maximal  element.   For  the  norm  of  the 
matrix  we  write 


(6.18)    ||M 


N 

/ 

> 

k=l 

J 

^J)^^  d^ 


=  M,   say. 
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With  this  notation  out  of  the  way  we  may  proceed  to  show 
the  contraction  mapping.   We  only  note  that  the  operator 
of  (6.l4)  maps  uniformly  continuous  and  bounded  functions 
into  the  same  space.   Further,  the  choice  of  an  initial 
iterate  is  relatively  unimportant  and  for  definiteness  we 
take  it  to  be  identically  zero.   Without  further  ado  we 
have , 


(6,19)  ||a"+1  -  a"  II  <  tMllA"  -  a"-1 


which,  for  t  <  T  <  |^  establishes  the  existence  and  unique- 
ness to  solutions  of  our  equation.   Further,  because  of  the 
norm  which  we  have  chosen,  the  iteration  converges  Itself  to 
a  function  which  is  uniformly  continuous.   Furthermore, 
since   M  is  independent  of   t   we  have  existence  and  unique- 
ness for  all  time.   For,  by  taking  the  solution  obtained  at 
t  =  T  as  the  initial  data  we  can  go  through  the  same  argument 
to  extend  the  solution  to  2T,  and  so  on.   We  can  use  this 
agrument  to  estimate  the  behavior  of  the  solution  in  time. 
We  choose   M   so  that  it  is  also  a  bound  for  the  initial  data 
given  by  (6.12),  then  from  (6.19) 


(6.20)         ||A^"^^^  -  a  II  <  MR  ,     for   0  <  t  <  T 


where 


(^•21)  R  =  i^ 
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On  passing  to  the  limit  In  (6.20) 


(6.21)         11^-^0  II  <  ^R  '     0  <  t  <  T 


Reapplying  this  argument  we  have,  in  general, 


(6.22)         ||A-A  II  <  Mr"      (n-l)T  <  t  <  nT 


or 


(6.23)         ll^~^o  I'  <  (M/R)R^/'^ 


so  that  the   A,   are  exponentially  bounded. 

So  far  our  remarks  have  been  restricted  to  the  moments 
A,  .   To  summarize,  we  have  shown  that  they  exist  for  all 
time  and  are  unique  in  the  class  of  uniformly  continuous 
and  bounded  functions,  and  that  they  are  exponentially  bounded 
in  time.   If  these  solutions  are  now  substituted  into  (6.9) 
it  furnishes  us  with  a  unique  solution  for  the  quantity  G. 
Further,  because  the  initial  conditions  in  (6.9)  have  the 
argument  (x-|t)  and  because  of  the  integral,   G  has  the  con- 
tinuous derivative  -S4r  +  i'    -s—   •   In  order  to  obtain  any 
further  information  on  the  differentiability  of  G   similar 
conditions  must  be  demanded  of  the  initial  conditions.   We 
have  been  able  to  obtain  all  of  these  results  under  the  mild 
boundedness  restrictions  imposed  on  the  initial  conditions  by 
the  remarks  in  connection  with  (6.12).   In  essence  only  the 
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boundedness  of  a  finite  number  of  velocity  space  moments 
were  required.   However,  as   N   Increases  (and  presumably 
equation  (6.1)  more  accurately  describes  the  Boltzmann  equa- 
tion), more  boundedness  conditions  are  required  of   G  . 

To  obtain  any  further  Information  on  solutions  of  model 
equations,  we  will  need  a  finer  analysis.   The  basis  for  this 
Is  furnished  by  the  dispersion  relation  analysis  which  was 
carried  out  earlier.   We,  therefore,  now  restrict  our  atten- 
tion to  one-dimensional  unsteady  equations.   Using  the  results 
of  the  previous  sections  we  can  show  the  uniform  boundedness 
of  solutions  and, moreover ,  display  the  asymptotic  character 
of  the  solution. 

In  order  to  accomplish  the  above  plan  we  must  further 
specify  the  model  equation  (6.1).   It  is  first  clear  from  the 
conservation  of  mass,  momentum,  and  energy  in  a  two  particle 
collision  that 

(6.25)  A^  =  A^  =  ^2  =  1 

o 
(when  t)^\,      1^1^  ^i>   t>2^^  ^    )•   Further,  from  the  construc- 
tion of  kinetic  models  given  in  reference   8  and  from  our 
normalization  the  remainder  of  the  A's  form  a  sequence  of 
non-negative  numbers  which  approach  zero. 

We  propose  to  solve  (6.1)  (actually  the  one -dimensional 
form  of  (6.1))  by  means  of  transforms.   Denoting  the  initial 
data  of  g   by  g  ,  the  transformed  system  can  be  written  as 


-  6^ 


N 
CO  ^   a  A  jz5 

where,  as  usual,  we  use  the  same  letter  to  denote  a  trans- 
formed quantity.   As  earlier  we  "reduce"  the  problem  by 
taking  moments  of  (6.26)  to  obtain 

N 

(6.27)         a  =  A  +  />   a  >^^C^^  ,     n  -  0,  ...,  N 
^    '  '  n    n   ^ — TT     m  m  mn 


m=0 


where 


(6.28)         A  = 


and 


(6.29)       C 

^        -" '  mn 


00 

031 


l+a-lk^. 


0  <  n,m  <  N 


00 


1+a-lk^, 
J  1 

—00 


C    Is  the  type  of  expression  found  In  Table  I.   An  Immediate 
mn  "^  -^       ^ 

property  of  C    which  will  prove  Important  later  is 


(6.30)         C   =  O(-)  ,    for  R  a  >  -1  and  \a\    -^' ^ 


This  is  true  for  any  bounded   k.   As  was  discussed  at  the  end 
of  Section  2,    Integrals  of  the  type  (6.28,  29)  define  two 
different  functions  depending  on  whether  k   is  greater  or 
less  than  zero.   We  shall  not  indicate  the  change  of  func- 
tional form  and  make  the  agreement  that  only  the  correct  form 
is  under  discussion.   The  change  of  functional  form  plays  no 
role  in  the  sequel  and  the  only  pertinent  remark  in  this  vein 
is  that  both  functional  forms  have  the  same  asymptotic  expan- 
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slon  (see  Appendix  l). 

The  linear  equation  given  by  equation  (6.27)  may  be 


Inverted  to  give   a  , 

n 
N 


(d-31)         a^  = 


y^  F    A 

^=^ — -  nm  m 


n    .   D 


where 


(6.32)         D  =  det(6    -AC   ),     0  <  m,n  <  N 


and  F    denotes  an  element  of  the  transposed  cofactor 
nm  ^ 

matrix  of   (6    -Ac   ).   We  shall  shortly  show  that  D  /^  0 
^  mn    m  mn^  ^  ~ 

for  5c  a  >  0   (which  at  this  stage  is  a  requirement  on  o) , 

permitting  us  to  perform  the  operation  of  (6.51)-   The 

vanishing  of  (6.32)  defines  the  dispersion  relation.   We 

shall  now  give  a  different  but  equivalent  representation  of 

the  dispersion  relation. 

Consider  the  following  relations: 

N 


GJ  2       b-^ 


(6.33)        q  =   "=^° 


n  n'  n 


1+a-ik^ 


(6.3^)   .       b^  =  /  j^^q  dl  n  =  0,  ..  .,  N. 


To  avoid  confusion  with  equation  (6.26)  we  have  replaced 

g  and   a   by  q  and   b  ,   respectively.   Upon  taking  the 

"  "  N      2 

moments  (6.34)  of  (6.33)  and  taking  ^  |b  |   /  0  we  obtain 

n=0   ^ 
an  eigenvector  b  ,   n  =  0,  1,  ...,  No   Equation  (6.32)  set 
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to  zero  defines  a  as  a  function  of  k  and  so  the  b  's 

n 

can  be  considered  as  functions  of  k  alone.   Next  sub- 
stituting these  expressions  for  b   Into  (6.35)  we  obtain 
q  as  a  function  of  k  and  i.      Alternately  using  this 
definition  of  q  we  automatically  satisfy  (6.3^).   Thus 
(6.33,3^)  are  equivalent  to  considering  a  as  a  function 
of   k  as  given  by  the  dispersion  relation  of  (6.32). 

We  now  take   q  as  having  been  obtained  by  the  above 
process,  and  further  we  take 


(6.35)         q  =  0)  'y~  b  j6 


Further,  letting 


(6.36)         Q  =  q/cD 


we  have  from  (6.33) 

(6.37)     a  /  CDQQd^  -  Ik  /  oje-^QQd^  =  -/  co(Q-  ^  b^A^j25^)Qde 


_N_ 

>_ 

n=0 


where  the  bar  denotes  the  complex  conjugate.   Taking  the 

real  part  of  (6.37) 

N 
-  /  a3(Q-  ^  h^T^JJq   d^ 
(6.38)         E  a   =   nzO_ 


n  n'^n' 


/  ojQQd^ 
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Prom  the  restrictions  placed  on  the  A   at  the  outset,  we 
have  (see  also  reference  'I  where  this  technique  was  applied 
to  get  similar  results) 

(6.39)         -1  <  9c  '^  <  0  • 


Prom  the  orthonormallty  of  the  6     we  have  the  condition 

J  ^n 


( 


2 


n= 
(6.^0)     <  or     2 


Q  =  2—  "^rfin  f  or  ^  a  =  0 

n=0 


Q  ^   ^  b^        for  5^  a  ;^  0 
n=0 


It  therefore  follows  from  (6.33)  that  if   a  =  0  then  k  =  0. 

We  now  wish  to  show  that  for   |k|  >  k   >  0,   that 
R.a  <  a  (k  )  <  0.   Otherwise  we  would  need  a  sequence  k. 
such  that,  as   k.  — >  k  ,   f^a(k.)  -^   0,  where  k   is  finite 
or  infinite,  which  is  a  contradiction  to  the  above  results. 

We  now  conclude  our  preparatory  work  by  some  crude 

estimates  on  D  and  P   .   From  (6.30)  and  under  the  same 

mn         \      ^    J 

restrictions  we  have  immediately 


(6.^11)         0-^1  +  0(l/a) 


and 


(6.42)         P   ^  1  +  0(l/o) 
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In  the  latter, unity  may  be  absent. 

Returning  to  our  problem  as  given  by  (6.27-29)?  we 
rewrite  the  Initial  conditions  as 


.2 
(6.45)         g^  -  e"^  [6^(k,^)  +  G^(k,e)] 


such   that 

(6.44) 

e    ^      G      =   0 

5 

k 

> 

e 

^   So 

} 

k 

< 

£ 

o 

To  avoid  tedious  arguments  we  assume  that   G  ,  G   are 
uniformly  bounded  In   ^   and   k,   and  further  that   G 
vanishes  for  |k|  sufficiently  large.   These  assumptions, 
especially  the  latter,  are  far  too  strong;  however,  their 
assumption  eliminates  messy  estimates.   The  quantity 
In  (6.44)  Is  still  left  open. 
We  first  show  that  for 


(6.45)         g^  =  e-^  G^(e,k) 


the  solution  to  our  problem  Is  exponentially  damped  In  time 
The  solution  for  a   Is  gotten  by  Inverting  a  .   Rather 
than  considering  the  entire  series  (6.31)?  we  can  consider 
a  typical  term,  which  we  denote  by  a  tllda. 
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(6.^6) 


a  - 

n 


(27r)2i 


dk  ■ 


B 


/'  -at  -Ik  „   /^  1  \ 
'  e    e    F  [o  ,k) 


mrj 


D(a,k; 


1+a-ik^ 


The   B  denotes  a  Bromwich  path  in  the  right  half  of  the 
complex  a-plane.   For   e  >  0   sufficiently  small,  we  can 
by  an  earlier  argument,  choose  a  6  >  0   such  that   D  will 
not  vanish  on   R-a  =  -6.   We  distort  our  original  path  B 
into  this  path,  ,^  of  Figure  6.1.   Further,  from  the  asymp- 
totic cases  studied  in  Sections  4  and  5>  we  know  that  we  can 
take 


(6.47) 


.2/2 


Moreover,  it  can  be  shown  (see  reference  h)    that  (6. ■4-7)  is 
sufficient  for  any  model  of  the  type  (6.1).   From  (6. ■4-1, 42^ 


we  may  choose  a   a    on 


■^^  such  that 


{  ( 

(6.48)    \a    I  <  e~^^  dk 
^         '  n 


mn 


D  •  a 


1^  '  1 '   m 
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G^(k,^)|d| 


J 


+ 


r 


( 


-6t 


I  dk   I 


-ia.t     -, 
e    ^  [1+0(^)1 

(l+a)[l+0(|)] 


e"^  /^G  (k,a)de 


1 


.k 


1+a 


-,  ,  a.  >a 
1  '  1 '   m 


where  we  have  written 


(6.49) 


a  =  -6  +  la. 
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The  Integration  of  the  first  term,  being  finite,  can  be 
majorized  by  a  constant   M-,   say.   In  the  second  term  the 
4  and  k  integrations  when  carried  out  are  seen  to  lead  to  a 
bounded  function  of  a.   Hence,  by  Dirichlet ' s  test  the 
integration  converges  and  is  majorized  by  Mp,  say.   We  have 
then 


(6.50)     |a^|  <  e"^^(M-^  +  M2)  =  e  ^^M 


It  is  important  to  note  that   M  is  at  worst  an  inverse 
power  of  £.    The  value  of  6  is  still  open  and  we  shall  fix 
it  in  the  next  phase  of  the  argument . 

We  now  go  through  the  same  type  of  analysis  for  initial 
conditions 


(6.51)     So  =  ^    G^(k,^) 


as  given  by  {6A^).      Again  It  is  only  necessary  to  consider  a 
typical  term.   Under  the  present  initial  conditions 


(6.52). 


n 


{2irfl 


dk 


00  f.  ^  , 
c   ^-at-lkx    ,^^^)  r  ^A   /2^  ,^^^)^^ 


J 


mn 


J 


D(a,k) 


1+a-ik^ 


-e   B 


As  before,   B  denotes  a  path  in  the  right  half  of  the  complex 
a-plane.   Since   e,   which  is  still  to  be  chosen,  is  arbi- 
trarily small  we  can  avail  ourselves  of  the  asymptotics 
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developed  In  the  earlier  sections.   We  first  exhibit  the 
asymptotic  expansion  of  the  last  Integral  of  (6.52).   Using 
the  identity 


(6.53) 


1       n  2 

=  1  +  X  +  X   + 


1-x 


n-1  ,  X 
+  n    + 


n 


1-x 


we  can  write 


(6.54) 


e 

1-  i^  i 

^        1+0     ^ 


d^  = 


^   «£^^  +  ITF 


4^/2- 


G  ^  dl  + 
e 


+ 


^  ^1+G^ 


n 


00    f,  d  / 
(     o-^  /2?n 


^^'G^d^ 


J 


ik 
1+a 


a.,  k 
°      1 
=  a  +  T^  + 
o   1+a 


k 


n 


+  o[ (itf)  1 


which  by  the  usual  arguments  shows  that  (6.5^)  is  an  asymp- 
totic development. 

In  considering  the  three  relaxation  model  we  showed 
from  the  asymptotic  analysis  that  for  k  — >  0,   there  are 
three  roots  to  D   in  the  neighborhood  of  the  origin  in  the 
complex  a-plane.   Each  of  these  roots  is  of  the  form  (see 
(5.18,  19)), 


(6.55)    a  =  iak  -  Sk^  +  O(k^)   , 


>  0 
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where  a   and  P   are  real.   It  is  furthermore  shown  In 
reference  4,  that  any  more  detailed  model  (given  in  general 
by  (6.1))  has  three  roots  in  the  neighborhood  of  the  origin. 
Since  the  values  of  k   are  restricted  by  the  integration  in 
(6.52)  we  may  enclose  the  three  roots  in  a  small  circle  about 
the  origin.   Further,  we  distort  the  path  B   so  that  it  is 
now  composed  of  the  circle  Xo   ^-^"^  the  imaginary  path  ^ . 
in  the  left  half  plane,  as  is  indicated  in  Figure  6.2.   Using 
the  estimates  used  previously  we  can  easily  show 


(6.56)   a^  = 


n 


27ri 


f 


dk 


-at-ikx -, 
e       da 


F  (a,k) 
m^  '  ' 


J 


D(a,k)(l+aj  ^ ^o^  1+a 


a-,  k 
[  a  +  ^  + 


-e 


C 


+  0(- 


k 


n 


{i+oY 


-)]    +   e"^^M 


where  the  latter  term  represents  the  contribution  of  the 

path  ^  -,  .   The  contribution  from  the  path  5-0   can  be 

evaluated  by  means  of  residues  to  given  a  contribution 

from  each  of  the  roots  of  D.   Using  the  representation  of 

a  typical  root  (6.55),  we  get  a  sum  of  three  terms  each  of 

which  may  be  written  as 

e 
f 

(6.57)    I  =    < 

-£ 

_  I 

It  is,  in  fact,  shown  in  the  cited  reference  that  any  finer 
model  will  not  alter  the  location  of  the  roots  to  O(k^). 


-Ikx+iakt-pk2t+0(k^)t^   +A,k+A^k2+...+0(k")]dk 

'-  o   1     2         V   /  J 
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where  A  's  are  constants  which  can  be  found.    Further,  by 
expanding  the  0(k  )  term  of  the  exponential  we  have 


(6.58)     I  = 


■lk(x-at)-pk  t^ (i+o(k5)t+. . . )(A^+. . .0(k")]dk 


-e 


We  now  choose 


(6.59)     e  =  O(i^)  , 


1  /    ^  1 
3  <  a  <  2 


With  this  choice  of   e   we  have  from  (6.47)  that  the  term 
e   M  which  we  obtained  In  connection  with  (6.46)  and  (6.52) 
Is  exponentially  decaying  In  time.   Estimating  the  remainder 
terms  we  get 


(6.60; 


and 


(6.61) 


e-P^  ^|0(k5)|t  dk  <  ^-^ 
'  ^   ^  '     ^  ,  3a-l 


e  ^    dk 


-e 


-Pk  t|^/,ni,|,,  ^ 
I 0(k  ) I dk  < 


na 


e  '^    dk 


-e 


where  the   c   Is  a  suitable  constant.   An  expansion  of  the 

0(k  )  term  shows  that  the  coefficient  of  the  Integral  In  (6.60) 


It  Is  possible  for  a   k  to  occur  In  the  denominator  of 
(6.57);  this  Is  easily  taken  care  of  by  a  principal  parts 
Integration. 
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-1/2 
can  be  taken  as  t  ^  .   To  complete  the  argument  which  shows 

that   I   of  (6.57)  Is  an  asymptotic  expansion,  we  observe. 


e    2 

r      -Pk  t  ^  n  ,, 
J  e  ^    k  dk 


(6.62)   ^:| 2 -^0(47o) 


■w^ 


-£ 


Finally  we  see  that  extandlng  the  Integration  to  Infinite 
limits  In  (6.58)  only  adds  an  asymptotically  negligible 
quantity  {^^    exponentially  decaying). 

For  the  purposes  of  this  study  we  have.  In  essence, 
proved  our  assertion  that  the  solutions  to  model  equations 
are  asymptotically  hydrodynamical .   In  Section  5  we  showed 
that  the  dispersion  relation  of  the  three  relaxation  time 
model  asymptotically  became  the  hydrodynamical  dispersion 
relation.   More  generally  It  can  be  shown  (see  reference  4) 
that  for  k  ^-'  0  and  o      In  the  neighborhood  of  the  origin, 
(6.33)  can  be  reduced  to  the  transform  of  the  Navler-Stokes 
equations  or,  more  generally,  to  some  higher  approximation 
of  the  Chapman-Enskog  procedure.   This  alone  is  not  enough  to 
give  us  the  assertion  stated  at  the  beginning  of  the  paragraph. 
This  is  so  on  two  counts.   First  (6.33)  involves  no  initial 
data  arid  secondly  we  have  no  license  to  Invert  this  system. 
The  above  analysis  removes  these  difficulties.   First  we 
have  shown  that  the  asymptotic  solution  only  involves   k —  0, 
0^0     which  by  the  previous  discussion  gives  the  Chapman- 
Enskog  terms.   Secondly,  the  initial  data  is  taken  care  of  by 
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(6.5^)  which  uses  only  moments  of  the  Initial  distribution 
function.    There  is  one  last  point.   Hydrodynamics, 
although  it  is  given  by  an  expansion  in   k  '^  0   in  kinetic 
theory,  regards  itself  as  exact.   It  then  requires  an 
integration  over  all   k   in  the  inversion.   However,  by  the 
last  paragraph  this  offers  no  difficulty. 

Integrals  of  the  type  (6.58)  lead  to  diffusing  waves 
in  general.   We  do  not  consider  them  further  in  this  study. 


Any  hydrodynamical  theory  only  uses  the  initial  data  of 
density,  temperature,  and  velocity.   It  is  interesting  to 
note  from  (6.5^)  that,  depending  on  the  asymptotic  order  of 
a  solution,  any  number  of  moments  of  the  initial  data  will 
enter.   This  will  be  taken  up  at  another  time. 

One  can  find  this  in  reference  10  where  the  fundamental 
solutions  of  the  linearized  Navier-Stokes  equations  are 
considered. 
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Appendix  I:   Evaluation  of  the  Fundamental  Integral. 

Each  element  of  the  matrix  C   of  Section  2  may  be 
expressed  in  terms  of  the  following  Integral  (see  Lemma  1 
of  Appendix  II) : 


(l.l)     M  -  (1+a) 


1+a-lk?. 


=  A 


00       „ 


J 


A-i 


\-  -   M(A) 


where 


(1.2) 


A  = 


l+g 

k 


As  mentioned  in  Section  2,    the  functional  form  of  (l.l) 
depends  on  the  sign  of  the  real  part  of  A,  in  any  case  for 
Re  A  ;^  0  we  have 


(1.3) 


M'  = 


dM 
dA 


A-i^. 


oxi^ 


-l^l)' 


From  this  it  is  easy  to  show  that 


(I.^) 


d   ^e-^'/^M 


dA 


A 


)  -  -e 


-a2/2 


+ 


Denoting  by   0~   the  limit  as  we  approach  the  origin  with 
Re  A  >  0   or  Re  A  <  o,   we  may  integrate  (l.4)  to  get. 


(1-5) 


A 


e  ""^  '^   dp.  +  const 
0^ 
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the  path  of  Integral  lying  In  the  appropriate  half -plane 
In  order  to  evaluate  the  constant  of  equation  (I.5)  we 
return  to  (l.l)  from  which  we  get 


+ 


(1.6) 


m; 

A 


+  =  1  nr^ 

A=0- 


slnce  the  principal  part  of 


(     code,. 


vanishes.   Making  use  of  the  error  function. 


(1-7) 


we  may  write 


^(x)  =  /2/7r       e   /   dt  (error  function) 


0 


/.„+ 


(1.8) 


M+(A)  =  YTJ2   A  e^  /^  (1  -  ^(A//2))   Re  A  > 


0 


M  = 


M 


p 
(A)  =  -)/V2  A  e^  /^  (1  +  l(A//2))   Re  A  <  0 


Using  the  asymptotic  properties  of  the  error  function 

+ 
we  now  list  the  asymptotic  properties  of  M~   that  will  be 

of  Interest  to  us. 
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(1.9)     M-(A)  =  1   -   l^+l^-ll>5+  ... 


+    (-1)  -^^^ ^+  0(-^, 


A 


for   I arg  A|  <  ^  -  e,  where   e  >  0  Is  small. 

i  +  i        2 
(1. 10)     M-(A)  =  {-if    "2^/27Ae^/2^^    1+1^ 

A"^    A^ 


10-5  ,      ,     (    .  ^n-l  l-3-5---(2n-3)  ,  ^fl   ^ 


for  jp  +  e  <  I  arg  A|  <  5^  -  £.    The  expansion  for  small 

values  of   A   is  immediately  found  using  the  error  function 

.   16 
expansion 


(1. 11)     i(A//2)  =  /277  {^    -^+   ^2 ^  -  •  •  •  ^ 

^       21  -2.5   3J  •2-^-7 
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Appendix  II:   Properties  of  the  Dispersion  Relation. 

We  will  now  prove  several  general  properties  of  the 
dispersion  relation  of  Section  2.   These  will  be  proven  in 
a  series  of  lemmas  given  below.   Actually,  the  properties 
exhibited  below  are  quite  general  and  apply  to  the  general- 
izations of  equation  (2.1)  given  in  reference  5«   Several 
of  the  arguments  are  for  convenience  repeated  from  reference 

Lemma  1.   Each  element   C. .   of  the  matrix  C   of  Section  2 

ij 

can  be  written  as 


(Il-l) 


,  n 


k 


where   F   is  a  real  function  of  A   which  is  either  odd  or 
even,  in  which  case   n   is  either  zero  or  one. 

We  see  from  Section  2  that  any  matrix  element  may  be 
written  as 


(II-2) 


''  cug(^^)d2,. 
1+0 -Iki-. 


where   g(^-,  )   is  either  an  odd  or  even  polynomial.   It 
suffices,  therefore,  to  consider 


(II. 3) 


CJJ^"  ^^1    i  (  '^^l  '^^l 

"  k 


1+a-ikt 


lA  +  t 


Using  the  binomial  expansion  we  may  write 
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n-1  . ^       / . ^ \n 


(II. 4)  ^J  =  (e^  +  lA)"  -  g^5^-^  lA  ...  -(lA) 

Now  consider  the  first  term  of  (3)  when  (4)  Is  substituted; 
this  gives 


(II. 5) 


co(^,  +  1A)"-1  d^ 


'1 


■  2 


-'-1/2 
which,  since  ao  oC  e   '^  ,   Is  an  even  polynomial  In   A   If 

n   Is  odd,  and  1  times  an  odd  polynomial  In   A   if  n   Is 

even.   By  Induction  each  term  of  (5)  contributes  the  same 

type  of  term.   And  hence,  since  g(£,-,  )  Is  either  an  odd  or 

even  polynomial  In  ^-.  ,    we  have  proven  the  lemma. 

In  addition  we  know  that  If  g  Is  an  even  polynomial, 
F  Is  odd  and  vice  versa.  An  Immediate  consequence  of  this 
Is  that  the  trace  elements  of  C   are  real  functions. 

The  above  lemma  states  that  each  element  of  the  deter- 
minant (2.20)  has  either  real  or  Imaginary  parity. 

Lemma  2 .   A  determinant  for  which  the  parities  of  any  row 
may  be  gotten  by  multiplying  the  preceding  row  by   1  (or 
any  determinant  which  may  be  put  Into  this  form,  say,  on 
extracting  1  from  a  row)  has  Itself  a  single  parity.   That 
Is,  the  determinant  expanded  Is  a  real  function  times   1 
or   1 . 

The  lemma  Is  clearly  true  for  a  two  by  two  determinant. 
Expanding  an  n  x  n  determinant  (which  satisfies  the 
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hypothesis)  along  the  first  column  we  get   n(n-l)  (n-l) 
determinants.   The  first  determinant  of  the  expansion 
satisfies  the  hypothesis,  the  second  determinant  can  be 
given  the  same  form  by  extracting  an   1   from  the  first 
row.   The  third  determinant  can  be  given  the  same  form  by 
Interchanging  the  first  and  second  rows.   In  fact,  each 
determinant  of  the  expansion  can  be  put  Into  the  form  of 
the  first  determinant  by  Interchanging  rows  and  extracting 
an   1   when  appropriate.   Imposing  the  Inductive  hypothesis 
we  prove  the  lemma. 

From  the  form  of  equation  12. l6)  and  Lemma  1,  we  see 
that  D(k,A)  satisfies  the  hypothesis  of  the  last  lemma. 
Therefore,  we  have 

Theorem  1:    '  If  D(k,A)  =  0,   then 

(II. 6)         D(k,A*)  =  0. 

This,  since   A  =  — —  ,   states  that  roots  of  the  dispersion 
relation  occur  In  conjugate  pairs.   With  the  use  of  Lemma  2, 
still  more  properties  of  the  determinant  may  be  shown. 

Lemma  3-   D   Is  of  even  degree  In   k   and   A. 

The  only  elements  of  C   odd  In  k  and  A   are  those  of 
Imaginary  parity.  '  But   D   Is  a  real  function,  hence  the 
terms  of  Imaginary  parity  must  occur  to  an  even  power,  and 
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the  lemma  is  proven. 
Finally  we  prove 

Theorem  2:   If  D(-k,-A)  =  0,   then  D(k,A)  =  0. 

It  Is  to  be  remembered  that  positive  and  negative  k 
lead  to  different  functional  forms  for  M,  denoted  by  M 
and  M~.   From  Appendix  I  we  have 

(II. 7)         M"(-k,-A)  =  M"^(k,A). 

Therefore,  denoting  the  corresponding  dispersion  relation 
by  the  same  subscripts,   if  k  >  0 


(II. 8)  D^(k,A)  =  D(k,A;  M+(k,?x))  =  0 


=  D(-k,-A;  M+(-k,-A)) 

which  by  (T)  proves  the  theorem.   For  the  purpose  of  this 
paper  this  says  that  we  can  restrict  attention  to  k  >  0 
or  k  <  0. 

Again,  since  A  =  ^^^  ,   the  latter  theorem  states 
that  +  k  lead  to  the  same  value  of  o. 
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TABLE    I 


C         -   ^ 


'12 


i   [1-M] 


C 
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'15 
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l/kA^ 
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0 


'21 
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